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Abstract. A generalization of the matrix model idea to quantum gravity 
in three and higher dimensions is known as group field theory (GFT). In this 
paper we study generalized GFT models that can be used to describe 3D quan- 
tum gravity coupled to point particles. The generalization considered is that 
of replacing the group leading to pure quantum gravity by the twisted product 
of the group with its dual —the so-called Drinfeld double of the group. The 
Drinfcld double is a quantum group in that it is an algebra that is both non- 
commutative and non-cocommutative, and special care is needed to define 
group field theory for it. We show how this is done, and study the result- 
ing GFT models. Of special interest is a new topological model that is the 
"Ponzano-Regge" model for the Drinfeld double. However, as we show, this 
model does not describe point particles. Motivated by the GFT considera- 
tions, we consider a more general class of models that are defined using not 
GFT, but the so-called chain mail techniques. A general model of this class 
does not produce 3-manifold invariants, but has an interpretation in terms of 
point particle Feynman diagrams. 



Recently a question of coupling of quantum gravity to matter has been revisited 
m works PE]. Both papers deal with quantum gravity in 2+1 dimensions, and 
analyze how matter Feynman diagrams are modified when the quantum gravity 
effects are taken into account. Both work contain results that are intriguing. How- 
ever, certain conceptual issues were left unclear. Thus, work pQ has suggested that 
the effect of quantum gravity is in modifying the measure of integration Y[ v dx v in 
a Feynman amplitude given by: 



Here one integrates over the positions of vertices x v of a Feynman graph, G(x — y) 
is the particle propagator in coordinate representation, and (vv') denotes edges 
of the Feynman diagram. In this scheme quantum gravity is only responsible for 
a modification of the Y[ v d>%v measure, and the particle content as well as the 
Lagrangian must be specified independently. In particular, one is not forced to 
using any particular propagator. Thus, this proposal is that of "quantum gravity 
+ matter". Work [2], on the other hand, showed that particular quantum gravity 
amplitudes can be interpreted as particle Feynman diagrams, provided one makes 
a special choice of the particle propagator. Thus, in this interpretation "quantum 
gravity = matter". These two interpretations seem to be in conflict. This paper 
is devoted to an analysis of these and related issues. We will argue that both 
interpretations are correct. 



1. Introduction 
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The basic philosophy that is to be pursued in the present paper is as follows. 
Consider the momentum representation Feynman amplitude of some theory. The 
precise definition of the theory is unimportant for us at this stage. We shall restrict 
our attention to general aspects common to Feynman amplitudes of all theories. 
Thus, the amplitude is given by: 

Z = J Y[ dp {vv/) G(p {vv/) )Y[s(^2p {vv/) ). (1.2) 

(vv') v v' 

Here the integrals are taken over the momenta P( vv ') on each edge of the digram; 
in addition, one has a momentum conservation law for each vertex. We have to 
modify the Feynman amplitude (|1.2|) so that it describes point particles. Point 
particles in 2+1 dimensions are conical singularities in the metric. As we shall see, 
this has the effect that particle's momentum is group-valued. For point particles in 
Euclidean 3-dimensional space, which is what we consider in this work, the relevant 
group is SU(2). The proposal to be developed in this paper is that point particles 
and processes involving them can still be described by Feynman diagrams l|1.2f) . 
but the momentum p on each edge of the diagram should be taken group-valued. 
The rest of the paper is devoted to various tests of this proposal. 

Some point particle theories we consider are topological; we refer the reader to 
the main body of the paper for the precise definition of what this means. As we shall 
see, these topological particle theories will be produced by the GFT. Other theories 
we consider are non-topological; GFT is of no direct relevance here. However, 
the techniques developed in the GFT context (such as the chain mail, see below) 
will prove extremely useful for these theories as well. This is how the GFT is 
omnipresent throughout the paper. 

We start by reminding the reader how point particles in 3D are described. For 
simplicity, this paper deals only with the case of 3 + O-dimcnsional spacetime, that 
is, with 3D metrics of Euclidean signature, and the cosmological constant is set to 
zero. 

1.1. Point particles. In 3 spacetime dimensions presence of matter has much 
more drastic consequences than in 4D: the asymptotic structure of spacetime gets 
modified. This has to do with the fact that the vacuum Einstein equations are much 
stronger in 3D: they require that the curvature is constant (zero) everywhere. Then 
an asymptotically flat spacetime is flat everywhere. When matter is placed inside 
it introduces a conical defect that can be felt at infinity. The modification due to 
matter is in the leading order, not in the next to leading as in higher dimensions. 
This has the consequence that the total mass in spacetime is proportional to the 
deficit angle created at infinity. Because the angle deficit cannot increase 27r, the 
mass is bounded from above. All this is very unlike the case of higher dimensions, 
and makes the theory of matter look rather strange and unfamiliar. But such 
strange features are at the end of the day responsible for a complete solubility of 
the theory. 

Point particles give the simplest and most natural form of matter to consider in 
3D. A point particle creates a conical singularity at the point where it is placed. 
The angle deficit 29 at the tip of the cone is related to particle's mass as 4irMG = 
9. Point particles move along geodesies. Their worldliness are lines of conical 
singularities. 
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To make this description more explicit, let us consider how a space containing a 
line of conical singularity can be obtained. This is achieved by performing a rotation 
of space, and identifying points related by this rotation. A rotation is described by 
an element g(29) 6 SU(2), which we choose to parameterize as follows: 

g (0) = e ^- s = cos6 + in- a sin9. (1.3) 

Here ft : ft ■ ft = 1 is a unit vector. It describes the direction of the axis of rotation, 
and thus the direction of particle's motion. The angle 29 specifies the amount of 
rotation, and thus gives the mass of the particle. The reason for choosing 9 and not 
29 for parameterizing the particle's mass will become clear below. An equivalent 
description in more mathematical terms is to say that a particle of a given mass 9 
is described by a group element g G Cg, where Cg is a conjugacy class in SU(2). 
Non-trivial conjugacy classes Cg, 9 ^ in SU(2) are spheres S 2 and a point in Cg 
describes the direction of motion. Thus, the information about both the mass and 
the direction of motion is encoded in g. This is why this quantity plays the role 
of momentum of a point particle. Note that the group-valued nature of particle's 
momentum automatically incorporates the bound on its mass. Note also that not 
only the upper, but the lower bound is naturally imposed as well - the SU(2) valued 
momenta are those of non-negative mass particles. 

Once points of R 3 related by a rotation g are identified, we obtain a space with a 
line of conical singularity in the direction ft. The angle deficit is 29. The space has 
zero curvature everywhere except along particle's worldline. One can compute the 
holonomy of the spin connection around the worldline and verify that it is equal to 
g(9). It is for this reason that the quantity 9 equal to half the angle deficit gives a 
more convenient measure parameterization of the mass. Let us now consider two 
point particles, with group elements 51,32 describing them. The holonomy around 
the pair is the product of holonomies, and the mass parameter 9 of the pair is given 
by: 

cos 6* = — Tr(<7i<?2) = cos6*i cos 9 2 — n~i ■ ft? sin^i sin #2- (1-4) 

Note that it does matter now in which order the product of holonomies is taken, for 
<7i<72 7^ 3231 as we are dealing with elements of the group. But in determining the 
mass of the combined system this ambiguity is irrelevant, for the trace is taken. Let 
us look at i|1.4|) in more detail. In case n~i is parallel to n 2 , this formula reduces to: 
9 = #1+6*2. When particles move in opposite directions we have: 9 = 1 6*i — 6*2 1 . Thus, 
the angle 9 is more properly interpreted as the magnitude of particle's momentum, 
not as its kinetic energy. 

It is instructive to compare 11.4|l to the usual law of addition of vectors in R 3 . 
Thus, consider two particles with momenta pi,p2. Define: m — \p\, ft = p/\p\- Then 
we have: 

m 2 — m\ + m 2 + 2n{ ■ n 2 m,]m 2 . (1-5) 

This formula for addition of momenta is to be compared with i|1.4|) . Indeed, as 
9 = AirMG, l|1.4|) can be expanded in powers of G. This gives exactly l|1.5|l . 
modified by 0(G 2 ) corrections: 

M 2 = Ml + Ml + 2n x ■ n 2 MiM 2 + (1.6) 

4ir 2 C 2 

(A/ 4 - Ml - M 2 4 - 3M X 2 M| - 4rTi • n 2 M 1 M 2 {Ml + M 2 2 )) + 0(G 4 ). 
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Thus, the fact that the law 1)1.4)1 of addition of momenta has the correct limit as 
G — > gives additional support to our identification of g with particle's momentum. 

Having understood why particle's momentum becomes group valued, let us re- 
mind the reader some basic facts about the Ponzano-Regge model 3 of 3D quantum 
gravity, and how point particles can be coupled to it. 

1.2. Ponzano-Regge model in presence of conical singularities. In Ponzano- 
Regge model, one triangulates the 3-manifold in question ("spacetime"), and assigns 
a certain amplitude to each triangulation. This amplitude turns out to be trian- 
gulation independent and gives a topological invariant of the 3-manifold. Let us 
remind the reader how the Ponzano-Regge amplitude 1)1. 11|) can be obtained from 
the gravity action. When written in the first order formalism, zero cosmological 
constant 3D gravity becomes the so-called BF theory. The action is given by: 

5[e,a] = f Tr(eAf). (1.7) 
Jai 

Here e is a Lie algebra valued one form, and f is the curvature of a G-connection 
a on M. The path integral of this simple theory reduces to an integral over flat 
connections: 

Z[M] = J PePa e iS [ e ' a J = J Pa 5(f). (1.8) 

A discretized version of this last integral can be obtained as follows. Let us trian- 
gulate the manifold M in question. Let us associate with every edge of the dual 
triangulation the holonomy matrix for the connection a. This holonomy matrix is 
an element of the group G. The product of these holonomies around the dual face 
must give the trivial element (because the connection is flat). Thus, one can take 
the ^-function of a product of group elements around each dual face, multiply these 
(5-functions and then integrate over the group elements on dual edges. This is the 
discretized version of 1)1.8)1 . By decomposing the S- function on the group G — SU(2) 
into characters, it is easy to show that this procedure gives the following amplitude: 

j2n&M3e)nm (i.9) 

It turns out, however, that this amplitude is not yet triangulation independent. In 
order to make it such one has to introduce a certain formal prefactor. Thus, let us 
consider the following sum over irreducible representations of SU(2): 

rr 2 :=^dim 3 2 . (1.10) 

3 

Note that this sum diverges and the quantity 7/ defined by it is formal. The sum 
above is equal to the 8- function on the group evaluated at the identity element. 
To make the amplitude l|1.9fl triangulation independent one has to consider the 
following multiple of it: 

r? 2v Eii dim ^)n( 6 i)» ( Ln ) 

{je} e * 

where V is the number of vertices in the triangulation used to obtain ()1.9(l . The 
amplitude 1)1.11(1 is now (formally, because of a prefactor that is actually equal to 
zero) triangulation independent. 
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The origin of this formal prefactor is now well-understood, see 0] . It has to do 
with the fact that the action (|1.7|l in addition to the usual SU(2) gauge symme- 
try has a non-compact gauge-symmetry. Indeed, because of the Bianchi identity 
(IaF — the action is invariant under the following transformation: B — > B + <Ia4>- 
Integration over B in the path integral includes integration over these gauge di- 
rections and produces infinities that has to be cancelled by the zero prefactor in 
(|l.ll|l . A more systematic procedure for dealing with these divergences was devel- 
oped in 0J|S] and requires a certain gauge fixing procedure. We shall not consider 
such a gauge fixing here and will proceed at a formal level, just keeping track of 
the (formal) pre-factors rj. All amplitudes can be rendered well-defined by passing 
to the quantum group context, that is, by considering the quantum group SU g (2) 
instead of SU(2). When q is a root of unity the cut-off on the spin present in the 
representation theory of SU 9 (2) makes the sum that defines r/~ 2 finite, and all the 
invariants well-defined. One should always have this passage to SU 9 (2) in mind 
when considering the ill-defined Ponzano-Regge model expressions. 

We have explained how the amplitude in (|1 - 1 1|> can be obtained from an integral 
over flat connections. It is now easy to add a point particle. Consider the case when 
a point particle is present and (a segment of) its worldline coincides with one of the 
edges of the triangulation of M. It will then create a conical singularity along this 
edge. The product of holonomies around the face dual to this edge is not trivial 
anymore. Instead it will lie in a conjugacy class determined by the mass of the 
particle (by its deficit angle). Thus, for dual edges that encircle particles' worldline 
one should put a (5-function concentrated on the conjugacy class determined by the 
particle. This gives a modified amplitude, with presence of the point particle taken 
into account. Let us introduce a 5-function that is picked on a particular conjugacy 
class 9: 



Here <5 e (x) is the usual 8- function on the group picked at the identity element. 
The sum in the second formula is taken over all irreducible representations of the 
group SU(2), and Xj{d) are the characters. The special group element hg is a 
representative of the conjugacy class 9 : hg = diag(exp i9, exp — iff). Note that 
when 9 = 0, the above (5-function reduces to the usual (5-function picked at the 
identity. Let us now assume that a point particle is present along all the edges of 
the triangulation, and denote the angle deficit on the edge e by 9 e . The modified 
quantum gravity amplitude in the presence of particles is given by: 



If particles are present along all the edges, there is no need to introduce a formal 
r\ prefactor as the amplitude is not required to be triangulation independent. If 
the particle is only present along some of the edges, one gets what can be called 
an observable of Ponzano-Regge model, see [7]. In this case pre-factors of 77 for 
vertices away from the particle edges are necessary to ensure a partial triangulation 
independence, see |7j for more details. 

As it was explained in and 0, the amplitude (| 1 . 1 3|> is the response of the 
quantized gravitational field to particle's presence. In the setting described the 
particle is present as en external source in the gravitational path integral. The 



Sg(x) := / dgS e (gh g g 1 x) = ^2xj(he)Xj(%) 



(1.12) 




(1.13) 



{Je} e 
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natural question is if one can integrate over the particle degrees of freedom and 
thus have a quantum theory of particles coupled to quantum gravity. In other 
words, the question is how to "second quantize" the particle degrees of freedom. 
Some sketches on how this can be done were presented in both 1 and |2J. What 
is missing in these papers is a principle that would generate Feynman diagrams for 
both gravity and particle degrees of freedom. This paper is aimed at filling this 
gap. As we shall see, the Feynman amplitudes for some of the theories (namely the 
topological ones) can be generated by GFT. To understand how this is possible, 
let us remind the reader some basic facts about the usual field theory Feynman 
diagrams, and, in particular, about a certain duality transformation that can be 
performed on a diagram. Let us remark that for the most part this paper deals 
with vacuum (closed) Feynman diagrams only. A generalization to open diagrams 
that are essential for doing the scattering theory will be left for future work. 

1.3. Momentum representation and duality. In addition to the coordinate 
(|l.l(l and momentum (|1.2(l representations of Feynman diagrams there is another, 
the so-called dual representation. Let us remind the reader how the dual repre- 
sentation is obtained in the case of two spacetime dimensions. To get the dual 
formulation we need to endow Feynman graphs with an additional "fat" structure. 
Namely, each edge of the digram must be replaced by two lines. The fat structure 
is equivalent to an ordering of edges at each vertex. The lines of the fat graph are 
then connected at vertices as specified by ordering. This structure is also equivalent 
to specifying a set of 2-cells, or faces of the diagram. Once this additional structure 
is introduced, one can define the genus G of the diagram to be given by the Euler 
formula: 2 — 2G = V — E + F, where V, E, F are the numbers of vertices, edges 
and faces of the diagram correspondingly. Note that the same Feynman diagram 
can be given different fat structures and thus have different genus. For example, 
the diagram with 2 tri-valent vertices and 3 edges (the (9-graph) can correspond to 
both genus zero and genus one after the fat structure is specified. 

Having introduced the fat structure, we can solve all the momentum conservation 
constraints. It is simplest to do it in the case of G = Q, so let us specialize to this 
situation. Let us introduce the new momentum variables, one for every face of the 
diagram. Let us denote these by pf. Given a pair (//') of adjacent faces, there 
is a unique edge (vv') that is a part of the boundary of both of them. A relation 
between the original momentum variables p/ V v>) an d the new variables pf is then 
given by: 

P{w>) = Pf ~Pf- (1-14) 

Let us check that the number of new variables is the number of old variables minus 
the number of constrains. The number of new variable is equal to F — 1: the 
number of faces minus one. We have to subtract one because everything depends 
on the differences only and does not change if we shift all the variables by the 
same amount. We have F — 1 = E — {V — 1) — 2G = E — (V — 1), where we 
have used our assumption that G = 0. Thus, the number of new variables equals 
to the number of the original variables, minus the number of 5-functions, equal 
to V — 1 because one 5-function is always redundant. In case G ^ one has 
to introduce an extra variable for each independent non-contractible loop on the 
surface, whose number is 2G. Thus, having expressed the momentum variables via 
Pf using (jl . 14|) - we have automatically solved all the conservation constraints. The 
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Feynman amplitude becomes: 



Y[d Pf n Gipj-pf,) 



(1.15) 



/ (//') 



We shall refer to (|1.15|) as the dual representation of a 2D Feynman diagram. 
One can perform a similar duality transformation in higher dimensions. The case 
relevant for us here is that of 3D, so let us briefly analyze what happens. Let us 
consider a vacuum (closed) Feynman diagram and introduce an additional structure 
that specifies faces, as well as an ordering of these faces around edges. For example, 
let us consider Feynman graphs whose edges and vertices are those of a triangulation 
of spacetime. As we shall soon see, such Feynman diagrams are natural if one wants 
to interpret the state sum models of 3D quantum gravity in particle terms. The 
faces of such "triangulation" diagrams are then the usual triangles. As in 2D, let 
us assign a variable pf to each face. The original edge momentum variable is then 
expressed as an appropriate sum of the face variables, for all the faces that share 
the given edge: 



One sums the face variables p / weighted with sign e/ ; a precise convention is unim- 
portant for us at the moment. The number of the original variables is E — (V — 1), 
and due to the fact that the Euler characteristic of any 3-manifold is equal to zero, 
we get: E — (V — 1) = F — (T — 1), which means that one in addition has to 
impose T — 1 constraints for the new variables. Thus, in the original momentum 
representation we had E edge variables together with the momentum conservation 
constraints for each vertex. In the dual formulation one has a momentum variable 
for each dual edge, and one conservation constraint for each dual vertex. We had 
to switch to the dual lattice in order to solve the 5-function constraints. 

1.4. Point particle Feynman diagrams. Here we further develop our proposal 
of modifying Feynman diagrams by making the momentum group valued. Thus, let 
us consider some theory of point particles in 3D that generates Feynman diagrams. 
A form of the action is unimportant for us. Our starting point will be diagrams 
in the momentum representation. As we have already described, the most unusual 
feature of point particles in 3D is that their momentum is group-valued and thus 
non-commutative. Thus, point particles in 3D behave unlike the standard relativis- 
tic fields in Minkowski spacetime. In spite of this, one can still consider Feynman 
diagrams. Indeed, in the momentum formulation (|1.2I) one only needs to specify 
what the propagator is, what the integration measure dp is, and what replaces the 
momentum conservation <5-functions. All this objects exist naturally for the 3D 
point particles. The propagator G(p) is some function on the group SU(2), which 
we shall leave unspecified for now. The integration measure is the Haar measure 
on the group. The momentum conservation is more subtle. The conservation con- 
straint becomes the condition that the product of the group elements <?(„„') for all v' 
is equal to the identity group element. However, it now does matter in which order 
the group elements are multiplied. Thus, some additional structure is necessary to 
define the momentum conservation constraints and Feynman amplitude as a whole. 
We shall specify this additional structure below. Ignoring this issue for the moment 




(1.16) 



8 



KIRILL KRASNOV 



wc get: 





9{vv') 



) 



(1.17) 



Here g e = gi V v') € SU(2) is the group element that describes particle's momentum 
on edge e, G{g) is a propagator which is a function on the group SU(2), and 5(g) 
is the usual 6- function on the group picked at the identity element. 

Let us now deal with the issue of defining the momentum conservation con- 
straints. It is clear that some additional structure is necessary. We have already 
seen an extra structure being added to a Feynman diagram when we considered 
the dual formulation. In that case an extra structure was introduced to solve the 
momentum constraints. It is clear that exactly the same structure can be used 
to define these constraints in the case of group-valued momenta. Thus, let us in- 
troduce an extra structure of faces, as well as ordering of faces around each edge. 
This is equivalent to introducing the structure of a dual complex. For simplicity 
we shall assume that the Feynman diagram we started from is a triangulation of 
some manifold. The dual complex is then the dual triangulation. Let us introduce 
new momentum variables: one for each face of the triangulation. The original edge 
momentum is expressed as a product of the face variables: 



The order of the product here is important and is given by the ordering of the faces 
around the edge: one multiples the holonomies across faces to get the holonomy 
around the edge. The new face momentum variables gf solve all the vertex mo- 
mentum conservation constraints. The amplitude (j 1 . 1 T|> in the dual formulation 
becomes: 



where TC are the tetrahedral constraints that were shown to be necessary by count- 
ing the variables in the previous subsection. 

1.5. Tetrahedron constraints. It turns out that the tetrahedron constraints have 
a simple meaning. 1 The constraints imply that the face variables are not all inde- 
pendent. There is a gauge symmetry that acts at each tetrahedron. At this stage 
it is convenient to pass to the dual triangulation. Recall that each face is dual to 
an edge e of the dual triangulation, and each tetrahedron is dual to a vertex v. 
To write down the action of the gauge transformation acting at v, let us orient all 
the edges e incident at this vertex to point away from v. Then the action on g e is 
as follows: g e — > g v g e , where g v is the same for all edges e incident at v. Fixing 
this gauge symmetry is the same as imposing the tetrahedron constraints in (|1 . 1 Q|) . 
Below we shall see how these constraints can be dealt with. 

There is another possible interpretation that can be given to the constraints. Let 
us consider a geometric tetrahedron in R 3 . The group valued face variables g t that 
we have introduced have the meaning of an SU(2) transformation that has to be 
carried out when one goes from one tetrahedron of the triangulation to the other. 

^We thank Laurent Freidel for helping us with this interpretation. 



9e = Y[ 9f 

f--eef 



(1.18) 




(1.19) 
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This group element can be represented as a product of two group elements, one 
for each tet: gj = {g t f)~ 1 g t f , where t, t' are the two tetrahedra that share the face 
/. Consider now: (g^) -1 ^,. This group element describes the rotation of one face 
/ into another /', and thus carries information about the dihedral angle between 
the two faces. Thus, the group variables gj introduced to solve the momentum 
constraints carry information about the dihedral angles between the faces of the 
triangulation. In particular, the total angle 2n — 9 e on an edge e is just the sum of 
all the dihedral angles around this edge: 

2n-e e =J2s f r, (1.20) 

<//') 

where the sum is taken over the pairs (//') of faces that share the edge e. This 
relation should be thought of as contained in the relation (|1.18|) . Having understood 
the (partial) geometrical meaning of the face variables as encoding the dihedral 
angles between the faces, we can state the (partial) meaning of the tetrahedron 
constraints. Recall that 6 dihedral angles of a tetrahedron satisfy one relation. 
This relation is obtained as follows. Consider the unit normals rii,i — 1, ... ,4 to 
all 4 faces of the tetrahedron. Let us form a matrix Ay = rii ■ rij of products of the 
normals. It is clear that ni ■ rij — — cos(0,y), and so the entries of Aj are just the 
cosines of the dihedral angles. However, as the vectors ni are 4 vectors in R 3 , they 
must be linearly dependent: 

^A«1 = 0. (1.21) 

i 

The coefficients A; here are just the areas of the corresponding faces. Because the 
normals fii are linearly dependent the determinant of the matrix A^ is equal to 
zero, which is the sought constraint among the 6 dihedral angles 6ij. This relation 
between the dihedral angles should be thought of as contained in the tetrahedron 
constraints in (|1.19|) . 

Let us now consider some example theories to which the above strategy can be 
applied. 

1.6. Pure gravity from point particle Feynman diagrams. One special im- 
portant case to be considered arises when particle's propagator is equal to the 
(5-function of the momentum. Had we been dealing with a usual theory in R 3 
in which the position and momentum representations are related by the Fourier 
transform such a choice of the propagator would mean having G(x — y) = 1 in the 
position space. Thus, this theory is a trivial one as far as particles are concerned - 
there are no particles. However, as we shall see in a moment, in the point particle 
context it leads to a non-trivial and interesting amplitude. 

From the perspective of the Feynman amplitude 11.17f) it is clear that this case 
is a bit singular. Indeed, the edge (S-functions guarantee that all the holonomies 
around edges are trivial. But then the vertex momentum conservation 5-function is 
redundant, and makes the amplitude divergent. To deal with this divergence, let us 
consider a modified amplitude with no momentum conservation 5-functions present. 
An equivalent way of dealing with this problem is to introduce the (formal) quantity 
i] given by (jl . 1()|> . Now we have to divide our amplitude by f]~ 2V , where V is the 
number of vertices in the diagram. We can now consider the dual formulation in 
which all the constraints are solved. Let us keep the factors of r], for they are going 
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to be an important part of the answer that we are about to get. The amplitude in 
the dual formulation is given by: 



Now it is easy to verify that the expression under the integral is invariant under the 
gauge transformations that act at the vertices of the dual triangulation. Thus, the 
tetrahedron constraints are redundant and can be dropped. This does not introduce 
any divergences, as the extra integrations performed if one drops the constraints 
gives the volume of the gauge group to the power of the number of tetrahedra. 
We always assume that the Haar measure on the group is normalized so that the 
volume of the group is one. Thus, when a gauge invariant quantity is integrated, 
the tetrahedron constraints can be dropped. 
Thus, let us consider: 



It is now possible to take all the group integrations explicitly Indeed, decomposing 
all the 5-function into characters, and performing the group integrations, one gets 
exactly the Ponzano-Regge amplitude with the important prefactor of rj 2V 

that is missing in l|1.9|) ■ It is only when this prefactor is present that the Ponzano- 
Regge amplitude is triangulation independent. 

Thus, we learn that the Ponzano-Regge amplitude (|1 . 1 1|1 . which, according to 
our previous discussion, is equal to the quantum gravity partition function, can 
also be interpreted as an amplitude for a point particle Feynman diagram with the 
particle's propagator given by the (^-function of the momentum. A few comments 
are in order. First, the Ponzano-Regge amplitude is independent of the triangu- 
lation chosen. Therefore, if we now take the particle interpretation, the Feynman 
amplitude is independent of a Feynman diagram. This could have been expected, 
because particle's propagator given by the (5-function essentially says that no par- 
ticle is present at all. This is why any Feynman diagram gives the same result. 
The second comment is on why such two different and seemingly contradictory in- 
terpretations are possible. Indeed, in one of the pictures, we are dealing with the 
pure gravity partition function, and no particles is ever mentioned. In the second 
picture, we consider point particles, whose momentum became non-commutative 
group-valued. This happened due to the back-reaction of the particle on the ge- 
ometry. Point particle Feynman diagrams as integrals over the group manifold 
take this back-reaction into account. Because pure gravity does not have its own 
degrees of freedom, the only variables we have to integrate over are those of the 
particles. Thus, particle's Feynman diagrams do give particle's quantum ampli- 
tudes with quantum gravity effects taken into account. This is why no separate 
path integration over gravity variables is necessary. 

The picture emerging is very appealing. Indeed, one has two different interpre- 
tations possible. In one of them one considers pure gravity with its "topological" 
degrees of freedom. The other picture gives a "materialistic" interpretation of the 
pure gravity result in terms of point particles. 




(1.22) 




(1.23) 
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It is now clear that a generalization is possible. Indeed, one does not have to 
take the propagator to be the <5-function. For example, one can consider the 5- 
function Sg picked on some non-trivial conjugacy class 9. The conjugacy classes 
can be chosen to be independent for different edges. Then in the dual formulation, 
one arrives at the Ponzano-Regge amplitude (|1.13(l modified by the presence of 
particles. 

The next level of generalization would be to allow all possible values of 9 on 
every edge, and integrate over each 9 e . This would mean integrating over possible 
masses of point particles on every edge. Unlike the theory with a given fixed set of 
angle deficits, which gives quantum gravity amplitudes modified by the presence of 
a fixed configuration of point particles, theory that integrates over particle's masses 
describes quantized particles together with quantized gravity. It is our main goal in 
this paper to study what types of theories of this type are possible. Thus, having 
understood that quantum gravity is described by taking into account a modification 
of particle's momentum that is due to the back-reaction, the question to ask is what 
fixes the theory that describes the particles. In other words, what is a principle 
that selects particle's propagator? What fixes valency of the vertices, that is the 
type of Feynman diagrams that are allowed? In the usual field theory we have 
the action principle that provides us with an answer to all these questions. In our 
case there is no such principle known. The models that we have considered so far 
are only known in the momentum representation or the corresponding dual. In 
fact, because the momentum space is non-commutative group manifold, one should 
expect the same to be true about the coordinate space. Thus, one can expect that 
the original coordinate representation of theories of the type we consider is in some 
non-commutative space. It has been argued recently that even in the no-gravity 
limit G — > this non-commutativity survives and leads to doubly special relativity, 
see |2j as well as a more recent paper and references therein for more details. 

Thus, one possible direction would be to understand what is the coordinate rep- 
resentation for point particle theories, and then postulate some action principle 
that would fix the propagator and interaction type. However, in this paper we shall 
instead consider the dual momentum formulation as fundamental. The question we 
address is which particle theories lead to "natural" dual representations. Indeed, 
to arrive to the quantum gravity interpretation in terms of the Ponzano-Regge am- 
plitude we had to switch to the dual representation. It is the dual representation 
that introduces a triangulation, and thus gives some relation to geometry. Another 
argument in favor of the dual representation is that Feynman diagrams in the orig- 
inal momentum representation are insensitive to the topology of the manifold they 
are on. Thus, they do not contain all of the degrees of freedom of the gravitational 
field, only the local point particle ones. To describe the global degrees of freedom 
one needs to have the extra structure that is available in the dual formulation. Fi- 
nally, even to define what the Feynman amplitude is in the case of a group-valued 
momentum requires the dual representation. Indeed, to define the momentum con- 
servation constraints one needs the structure of faces of the diagram, which is only 
available with the dual graph. All this makes it rather clear that the dual formu- 
lation is indispensable for the description of point particles in 3 dimensions. The 
logic of this paper will be to start from the formulation in terms of group field 
theory and the dual graph and then give the Feynman diagram interpretation. 
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With this idea in mind, we have to discuss how to obtain theories for which a 
dual formulation is available. A natural way to do this is using a generalization of 
the matrix model idea. Such a generalization is available already for the Ponzano- 
Regge model, and is known under the name of Boulatov theory We shall 
describe similar theories that incorporate point particles. Before we do this, let us 
remind the reader some basic facts about the usual matrix models. 

1.7. Matrix models. Much of the activity in theoretical physics at the end of 
1980's beginning of 1990's was concentrated in the area of matrix models. The idea 
is that the perturbative expansion of a simple matrix model given by the action: 



where M is, say, a hermitian N x N matrix, and V(M) is some (polynomial) 
potential, can be interpreted as a sum over discretized random surfaces. Such a 
sum, in turn, gives the path integral of Euclidean 2D gravity. The above matrix 
model can be solved by a variety of techniques and the predicted critical exponents 
match those obtained by the continuous methods. The simple one matrix model 
given above corresponds to the theory of pure 2D gravity, or, if one chooses the 
potential V(M) appropriately, gives 2D gravity coupled to the so called (2, 2m — 1) 
matter. Couplings to other types of c < 1 matter are possible by considering the 
multi-matrix models, see e.g. ||] and references therein for more details. 

More generally, one can consider a matrix model for an infinite set of matrices 
parameterized by a continuous parameter that is customarily referred to as time t. 
The corresponding matrix quantum mechanics gives 2D gravity coupled to a single 
scalar field, i.e. the so-called c = 1 matter, and is given by the following action: 



Here the matrices M(t) are functions of time, and to is a "mass" parameter. The 
matrix quantum mechanics (|1.25(l can be solved exactly. The results can be com- 
pared with those obtained in the continuous approaches, with full agreement, see 
[5j for more details. 

1.8. Field theory over a group manifold. The idea of generating random 
discretized manifolds from a matrix model was generalized to 3 dimensions in 
PH- Here one considers a field 4>{gi, 92,93) on three copies of a group mani- 
fold G. The field 4>{gi, 92, 93) should in addition satisfy an invariance property: 
4>{9i9, 929, 339) = 4>{gii 92, 93)- One should think of each argument of the field as a 
generalization of the matrix index in 2D. One considers the following simple action: 



Here the first product in the interaction term is over pairs of indices i, j = 1, . . . , 4, 
and indices i,j, k, I in the arguments of <j> are such that i =/= j k I. Using the 
Fourier analysis on the group one can expand function </> into the group matrix 
elements. One then obtains the Feynman rules from the action written in the 
"momentum" representation. The resulting sum is over triangulated 3-manifolds, 
or, more precisely, pseudo-manifolds. Edges of the Feynman graph correspond 
to faces of the triangulation, and Feynman graph vertices corresponding to the 




(1.24) 




(1.25) 




(1.26) 
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tctrahcdra. Edges of the complex get labeled by irreducible G-representations l e 
and one has to sum over these. The perturbative expansion is given by: 

e .g T En dim ^)ii^)- ( l27 ) 

complexes Oe} e * 

The first product here is over the edges of the triangulation, and the second is over 
its 3-cells (tetrahedra) . The quantity T is the total number of tetrahedra in a com- 
plex. The quantity (6j) is the Racah coefficient that depends on 6 representations. 
The amplitude of each complex here is (almost, because of the missing rj factors) 
the Ponzano-Regge amplitude |3] . We shall give details of the calculation that leads 
to this result in the next section. 

Thus, the generalized matrix model of Boulatov [Hj gives a sum over random 
triangulated 3-manifolds weighted by the exponent of the discrete version of the 
gravity action. Unlike the case of 2D matrix models it was not possible to solve 
the Boulatov theory, and not so much has been learned from it about 3D quantum 
gravity, see, however, for an interesting proposal for summing over the 3- 
topologies. Similar theories have been considered in higher number of dimensions 
as well. This more general class of models defined on the group manifold and 
leading to sums over triangulated manifolds have received the name of group field 
theories (GFT). 

1.9. Outline and organization of the paper. We have seen that the Ponzano- 
Regge amplitude has a point particle interpretation: it is the amplitude for a theory 
of point particles with a propagator being the J-function in the momentum space, 
and the amplitude is computed in the dual momentum representation. In this 
paper we shall address the following question: what kind of group field theory can 
give dual momentum representation for point particle theory with a non-trivial 
propagator? As we shall see, the requirement that Feynman diagrams generated 
by group field theory also admit an interpretation as dual of some other Feynman 
diagrams severely restricts the possible choices of theories. There is essentially a 
few possible choices, which we shall discuss in due course. 

The main idea of this paper is to consider a group field theory on the Drin- 
feld double D(SU(2)) of the group SU(2). To the best of our knowledge, the fact 
that point particles in 3D are naturally described by D(SU(2)) was discovered and 
explored in More recently, it was used in to give the CS formulation of 

the Ponzano-Regge model. As we shall see, a group field theory construction on 
D(SU(2)) similar to that of Boulatov [Hj can be performed. In some sense, the idea 
of going from SU(2) group field theory to D(SU(2)) one is analogous to the way 
one introduces the matter degrees of freedom in 2D, where multi-matrix instead 
of one-matrix models are considered. However, the analogy here is not direct, for, 
as we shall find, the most obvious generalization does not lead to a theory with 
point particles in it. Nevertheless, group field theory considerations will allow us to 
develop certain methods that will later be used to describe a general point particle 
theory. 

The organization of this paper is as follows. In the next section we describe 
how Ponzano-Regge model is obtained from Boulatov theory. Here we also give an 
abstract algebraic formulation of Boulatov model that will be later used to define 
GFT for the Drinfeld double. We describe the Drinfeld double of a classical group 
in sectional Here we discuss irreducible representations of D(SU(2)), characters, as 
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well as an important projector property satisfied by them. Group field theory for the 
Drinfcld double is defined in section^] and some of its properties are proved. Section 
defines more general particle models that are not related to GFT. Interpretation 
of all the models is developed in section We conclude with a discussion of the 
results obtained. 

Let us note that, apart from the works already mentioned, the question of cou- 
pling of quantum gravity to matter was considered in paper |12) . However, the 
approach taken in the present paper is rather different. 



2. BOULATOV FIELD THEORY OVER THE GROUP 

Boulatov theory [H] can be formulated in several equivalent ways. One of this 
formulation will be used for generalization to the Drinfeld double. Let us first give 
the original formulation used in [5] 

2.1. Original Boulatov formulation. Let us consider a scalar field on the group 
manifold </>(g). For definitcness we take the group to be G — SU(2) that corresponds 
to 3D Euclidean gravity. However, one can take any (compact) Lie group. The 
resulting theory would still be topological, and be related to the so-called BF theory 
for the group G. The Fourier decomposition of the function (j) is given by: 

<Ks)=EE*m«(<?)Cn- (2- 1 ) 
; ran 

Here t mn (g) are the matrix elements in Z'th representation, and <j) mn are the Fourier 
coefficients. It is extremely helpful to introduce a graphic notation for this formula. 
Denoting the Fourier coefficients by a box, and the matrix elements by a circle, 
both with two lines for indices m, n sticking out, we get: 

m = E ^ 1 (g ) (2.2) 

The sum over m, n is implied here. The graphical notation is much easier to read 
than 12. and we shall write many formulas using it in what follows. 

The field of Boulatov theory is on 3 copies of the group manifold. Analogous 
Fourier expansion is given by: 



Hgi,92,g 3 )= £ r^l YgS (t 2 ) (t) (2.3) 



The field is required to be symmetric: 

0(9<r(i), 9a(2),9a(3)) = (~ 1) <K#i ) 92 , 53) ■ (2.4) 

Here a is a permutation, and \a\ its signature. In addition, it is required to satisfy 
the following invariance property: 

<t>(9i9, 929, 535) = <j>(gi,92, 93)- (2.5) 

Let us find consequences of this for the Fourier decomposition. Let us integrate the 
right hand side of (|2.5(l over the group. We are using the normalized Haar measure, 
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so we will get 4>(gi, 52, 53) on the right hand side. On the left hand side we can use 
the formula H7.3(l for the integral of three matrix elements. Thus, we get: 



^(31,52,53)= Y 




(2.6) 



Let us now introduce a new set of Fourier coefficients <t> l m?m 2 m 3 ■ Graphically, they 
are defined by: 



(2.7) 



We shall only use the new, modified set of coefficients <|2.7|) till the end of this 
section. Thus, we shall omit the tilde. Our final expression for the Fourier decom- 
position is given by: 



0(51,52,53) 




(2.8) 



It is now straightforward to write an expression for the action in terms of the 
Fourier coefficients. The action (|1.26fl is designed in such a way that there is always 
two matrix elements containing the same argument. Thus, a repeated usage of the 
formula (|7.2|l gives: 



SI 



£1/2^3 



h_ 
h 



(2.9) 
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-M34 




<P 





Here the first sum is taken only over sets h,h,h that satisfy the triangular in- 
equalities, and this is reflected in a prime next to the sum symbol. To arrive to 
this expression we have used the normalization of the intertwiner given by l|7.4|l . 
We have also introduced the so-called Racah coefficient, or the 6j-symbol that is 
denoted by two rows of spins in brackets. 

It is now easy to derive the Feynman rules. The propagator and the vertex are 
given by, correspondingly: 



I 



(2.10) 
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where the black box denotes the sum over permutations, and 

hi/ 

( I12 ha hi \ — 

\ ^34 ^24 hs / 124 

123 \ 




(2.11) 



Note now that each Feynman graph of this theory is a dual skeleton of some sim- 
plicial 3D manifold. Indeed, vertices of the graph correspond to simplices, edges 
corresponds to faces, and faces (closed loops) corresponds to edges. A proper way 
to describe the 3-manifolds arising is in terms of pseudo-manifolds. It will be pre- 
sented below. Each Feynman amplitude is weighted by almost the Ponzano-Regge 
amplitude. What is missing is a (vanishing for the classical group) factor of rj 2 for 
each vertex of the triangulation. 

Note that, unlike the case of 2D gravity where we have a clear interpretation 
of both the rank TV of the matrices and the coupling constant g (they are related 
to Newton and cosmological constant correspondingly), in the case of Boulatov 
theory the interpretation of g is obscure. It is not anymore related to the volume 
of a tetrahedron, for the latter is now a function of the spins. The appearance of 
Newton's constant is also not that direct. It serves to relate the spin l e labeling 
edges to their physical (dimensionful) length. It has been argued recently ^3] that 
the coupling constant g should be thought of as the loop counting parameter that 
weights the topology changing processes. However, this issue is not settled and we 
shall not comment on it any further. 

2.2. Algebra structure. Let us remind the reader that the algebra A* of functions 
on the group can be given a structure of the Hopf algebra. The reason why the 
algebra is referred to as A* and not A will become clear below. 2 The Hopf algebra 
structure is as follows: 

multiplication (/i * f 2 ){x) = fi(x)f 2 (x), 

identity i(x) = 1, 

co — multiplication (A*f)(x, y) = f(xy), 

co -unit e(f) = f(e), (2.12) 

antipode (nf){x) = f(x~ l ), 

involution f°{x) — f[x). 

Of importance is also the Haar functional h : A — > C: 

h*(f) = J dxf(x). (2.13) 

A non-degenerate pairing 

(/i,/ 2 > = / dxh(x)h{x) (2.14) 
identifies A* with its dual (A*)* = A. 



2 tt is customary to denote the commutative point-wise product of functions by • and the non- 
commutative convolution product by *, so we go against the convention here. The reason for this 
unusual notation is that we decided to embed A into the Drinfeld double D(SU(2)) in a particular 
way, see below, and this induces the • and * products in the way chosen. 
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The dual algebra A is also a Hopf algebra. A multiplication on A is denoted by 
• and is introduced via: 

</i»/2,0>:=(/i®/a,A* 5 >. (2.15) 

One obtains: 

(jW 2 )(z) = J dzf 1 {z)h{z- 1 x). (2.16) 

All other operations on A read: 

identity i(x) = S e (x), (2.17) 
co - multiplication (Af)(x, y) = f(x)5 x (y), 

co -unit £(/) = J dxf(x), (2.18) 

antipode (nf)(x) = /(x _1 ), 
involution S*( x ) = /(.t _1 ). 
It is easy to check that * is indeed an involution: 

(A • S2T = Si • A*- (2.19) 

One also needs a Haar functional h : A — ► C: 

M/) = /(e)- (2.20) 
Using this structure, a positive definite inner product can be defined: 



(fuh) = J dxSx{x)f 2 {x) = h(Sx • St) = KSZ • Si)- (2-21) 

The algebraic structure on A will play an important role when we give an algebraic 
formulation of the Boulatov theory. 

2.3. Projectors. Of special importance are functions on the group that satisfy the 
projector property: 

/•/ = /• (2.22) 

Examples of such projectors are given by: (i) characters 

: d:m,\, ;•<•*• (2.23) 

and (ii) spherical functions dimjTQ (a;). In both cases the projector property is 
readily verified using l|7.1|l . We have: 

P j *P k = 8 jk P j . (2.24) 

We will mostly be interested in these character projectors. We note that the unit 
element, that is, the (^-function on the group can be decomposed into the character 
projectors: 

i(x)=S e (x)=J2 P i- (2-25) 

j 

Of importance for what follows is the following element of A* ® A* : 

G = (k ® id)A*i, (2.26) 
G(x,y) = 5 e (x~ 1 y). 
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In view of a property: 

' dyG(x,y)G(y,z) = G(x,z) (2.27) 



this function on G x G can be referred to as a propagator. The propagator is just 
a kernel of the operator G acting on A via the "-product. 

The unit element and the characters are the basic projectors on a single copy of 
A. More interesting projectors can be constructed when working with several copies 
of the algebra. Consider an element 9 G A® 3 obtained from three propagators G: 

0:= (l m ®h 3 )G<S)G<S)G. (2.28) 

Here we have introduced a new operation h-$ : A* ® A* ® A* ^ C: 



M/1,/2,/3) = &*(/i*/ 2 */s) = / dxh(x)f 2 (x)h(x). (2.29) 
As a function on 3 copies of the group the clement 9 is given by: 



where 



J dgS e (x 1 1 g)S e (x 2 1 g)S e (x 3 1 g) = ^J^dim^^^), (2.30) 



9 n (x l ) := e jld2tj3 (xi,x 2 ,x 3 ) := / dg x n ( x i 1 9)Xj 2 ( x 2 1 9)Xj 2 { x 2 V)- ( 2 - 31 ) 
Let us note that 9 is ^-invariant: 9* — 9. It is easy to verify that 9 is a projector: 

6 •6 = 6. (2.32) 

To show this one uses (17.31 and l|7.4|) . The function 6{x\, x 2 , £3) projects onto 
functions invariant under the left diagonal action. Indeed, consider an arbitrary 
function 4>{x\, x 2 , £3) G A® 3 . Define the <E> to be 9 with the projector 9 applied on 
the left. The function <!> is invariant under the left diagonal shifts: 

(j>:=6»4>, 4>(gxi,gx 2 ,gx 3 ) = 4>(xi,x 2 ,x 3 ). (2.33) 

The field <j> obtained this way is the basic field of Boulatov theory, see (|2.5|l . The 
above more abstract formulation in algebra terms in necessary for a generalization 
to field theory on the Drinfeld double. 

2.4. Algebraic formulation of Boulatov's theory. An equivalent formulation 
of Boulatov theory can be achieved using the Hopf algebra operations introduced 
earlier in this section. More specifically, we will need the second algebra structure 
in terms of the non-commutative ^-product and *-involution. Let us first write the 
kinetic term of the field theory action l|1.26fl . We have: 

Sun = \h{{9 • • (9 • 0)) = ±h(4>* •6*4>). (2.34) 

Here we have used the fact (I2.19f) that * is an involution, and the fact that 9 is 
a projector. All operations here, namely the ^-product, the ^-involution, and the 
co- unit are understood in this formula as acting on A® 3 , separately on each of the 
copies of the algebra. 

There are two possible points of view on the kinetic term Ij2.34|l . One is that the 
action is a functional of a field <fi that satisfies the invariance property 9 • 4> — 4>- 
This is the point of view taken in the original formulation of Boulatov theory that 
has been described above. The other interpretation is suggested by Ij2.34|l . and is to 
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view 9 in this expression as the kinetic term "differential operator" , or the inverse 
of the propagator of the theory. The action is then invariant under the following 
symmetry: S(f> = ip, where ip is a field satisfying 9»tp = 0. This symmetry should be 
viewed as gauge symmetry of the theory; the "physical" degrees of freedom are not 
those of 4> but those of the 0-cohomology classes. Boulatov theory in its original 
formulation H1.26[l can be viewed as the gauge-fixed version of this theory. The 
gauge symmetry described is of great importance. For example, one can consider 
a theory with no ^-projector inserted in the action. This theory does not have any 
gravitational interpretation, as it is easy to check. Thus, it is the presence of 9 in 
the action, and thus the extra gauge symmetry that ensures a relation to gravity. 
It can be argued that this symmetry is the usual diffeomorphism-invariance of a 
gravitational theory in disguise. 

Once the kinetic term is understood, it is easy to write down the interaction 
term as well. One has to take 4 fields 9 • <f> and •-multiply them all as is suggested 
by the structure of the interaction term in (|1.26(l to obtain an element of A® & . 
Each copy of the algebra is a product of two fields; the ^-involution has to taken 
on one of them. After all fields are multiplied, the Haar functional h has to be 
applied to get a number. We will not write the corresponding expression as it is 
rather cumbersome. The structure arising is best understood using the language of 
operator kernels to which we now turn. 

2.5. Formulation in terms of kernels. We have given an algebraic formulation 
of Boulatov theory, in which fields are viewed as operators in A® 3 , and one uses 
the "-product on A to write the action. An equivalent formulation can be given by 
introducing kernels of all the operators. Such a formulation is more familiar in a 
field theory context, and will be quite instrumental in dealing with the theory. 
Let us interpret the formula (|2.16(l as follows. The quantity: 

h{x,y) := f2(x- 1 y) = ( K ®id)A*/ 2 , (2.35) 

which is an element of A* ® A* , is interpreted as the kernel of an operator Of 2 
corresponding to that acts on j\ from the right: 

fi°O h = J dzfx(z)f 2 (z,x). (2.36) 

Thus, instead of dealing with operators from A one can work with kernels from 

A*® A*. 

The inner product Ij2.21|l can also be expressed in terms of the kernels. Since 
f(y, z) := f(y~ 1 z), the kernel for /* is given by: 

f*(y,z) = J{^). (2.37) 

Therefore: 

(fi,h) = J dxdy f 1 (x,y)f^(y,x). (2.38) 

Here we have taken a convolution of the kernels of two operators, and then took 
the Haar functional given by the trace of the corresponding kernel: 



Kf) ~> / dxf(x,x). 



(2.39) 
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Let us now introduce kernels for all the objects that are necessary to define 
Boulatov theory. The kernel for the field cj> 6 A* ® A* ® A* is: 

:= (( K ®id)A*)® 3 0, (2.40) 
<£(a:i, 3/152:2,2/2; 2:3,2/3) := ^(x^yi, x 2 ~ 1 y 2 , x^ya). 
The kernel for the projector is similarly given by: 

d(xi,yi;x2,y 2 ;x 3 ,y 3 ) := flfc^V, z^ 1 ^, 2:3 1 2/3)- (2-41) 
The action of 6 on the field now reads: 

$(2:1, 2/1; 2:2, 2/2; 2:3, 2/3) := (2.42) 
J dz 1 dz 2 dz 3 (j)(x 1 ,z 1 ;x2, z 2 - 1 x 3 ,z 3 )e(z 1 ,y 1 ; z 2 ,y 2 ; z 3 ,y 3 ). 

It is now easy to write the action for the theory. We shall assume that the field 
is real. The action reads: 

1 f 3 

S = - Y^dxidyidzi <f>(xi,yi)6(yi,Zi)(f>(zi,Xi) + (2.43) 

i— 1 

f 4 

tt / J~~J dxijdyijdzij cj){xij , y^j 

Here x^ = Xji , Zij — Zji , but j/y ^ t/^j and are independent integration vari- 
ables. We have also introduced a notation: <fi(xi,yi) := 4>(xi,yi; x 2 ,y 2 ;x 3 ,y 3 ) and 
(p(xij,yij) := (f>(xij,yij;Xik,yik;Xii,yu) with ijLj^k^l. 

2.6. Remarks on Boulatov theory. Let us remark on possible choices of prop- 
agators in 12.43|) . One can try to use a different operator in place of 9 in both the 
kinetic and the potential terms. Or, more generally, one can have an interaction 
term that is built of quantities different than the one used in the kinetic term, see 
|14| where such more general theories are considered. However, this more general 
class of theories fails to lead to 3-manifold invariants. Related to this is the fact that 
these more general theories do not have a dual Feynman diagram interpretation. 
As we have discussed in the introduction, we would like to restrict our attention to 
a special class of theories, namely those in which Feynman amplitudes that follow 
from the group field theory expansion also have an interpretation in terms of Feyn- 
man diagrams for the dual complex. In other words, as we have discussed, group 
field theory Feynman diagrams have the interpretation of amplitudes for a complex 
dual to some triangulated 3-manifold. The triangulation itself can be considered 
as a Feynman diagram. When can the group field theory amplitude, which is the 
one for the dual triangulation, be interpreted as a Feynman amplitude for the orig- 
inal triangulation as well? Thus, the question we are posing is which group field 
theories admit a dual formulation. As is clear from our discussion of Boulatov and 
Ponzano-Regge models in the introduction, this particular theory does admit both 
interpretations. Are there any other theories with a similar property? 

We shall not attempt to answer this question in its full generality. We shall 
only make a remark concerning the choice of the propagator of the model. As 
we have seen, the propagator of any 3d group field theory model should consist 
of 3 strands. Thus, it can be described as made out of 3 propagators, one for 
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each strand. This is what happens in the case of Boulatov model, where the 9- 
projector is made out of 3 G-propagators (|2.26(l . Now each strand of the group 
field theory diagram that forms a closed loop is interpreted as dual to an edge 
of a triangulated 3-manifold. When this triangulation is itself interpreted as a 
Feynman diagram, there must be a propagator for every edge. It is clear that this 
propagator will be built from the group field theory strand propagator G, and will 
just be a certain power of it, with the power given by a number of dual edges 
forming the boundary of the dual face. Because different triangulations can have 
this number different, to have the interpretation we are after the propagator G 
must be a projector G • G = G. Only such propagators admit the dual Feynman 
diagram interpretation. This requirement is clearly very restrictive, and limits the 
choice of possible propagators dramatically. In this sense the models of the type 
we are considering in this paper are very scarce. 

Now that we have understood how to formulate Boulatov's theory in abstract 
terms, let us construct an analogous theory with the group manifold replaced by 
the Drinfeld double of SU(2). 



The Drinfeld double of a classical group was first studied in ^S] ■ Our description 
of the Drinfeld double D(G) of G = SU(2) closely follows that in [Ej. Following 
this reference, we describe D(G) as the space of functions on it. As a linear space 
D(G) is identified with the space C(G xG) of functions on two copies of the group. 
On D(G) we have a non-degenerate pairing: 



This pairing identifies the dual D{G)* of the Drinfeld double with C(G x G). The 
following operations are defined on D(G): 



3. Drinfeld double of SU(2) 




(3.1) 



multiplication 




identity 
co — multiplication 



l(a;,3/) = S e (y), 

(A/Xxi, yv,x 2 , 3/2) = f(xix 2 ,yi)S yi (y 2 ), 



co — unit 




(3.2) 



By duality we have the 



involution 



antipode 



(Kf)(x,y) = f(y 1 x 1 y,y x ), 
f*(x,y) = f{y~ 1 xy,y- 1 ). 
following operations on the dual D(G)*: 



multiplication 




identity 
co — multiplication 



i(x,y) = S e (x), 

(A*f)(x 1 ,y 1 ;x 2 ,y 2 ) = f(x 1 ,y 1 y 2 )S X2 (y 1 ; 1 x 1 y 1 ), 



co — unit 




(3.3) 



antipode 
involution 



(K*f)(x,y) = f(y l x 1 y,y *), 
f°(x,y)=f(x-\y). 
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The universal R-matrix R 6 D{G) <8> D(G) is given by: 

R(xi,yv,x 2 ,y 2 ) = 5 e (yi)S e (x 1 y 2 ~ 1 ). (3.4) 
We will also need the central ribbon element: 

c(x,y) = S e (xy) (3.5) 

and the monodromy clement: 

Q(xi,yi;x 2 ,y 2 ) ■= (R»aR)(xi, yi;x 2 ,y 2 ) = 5 yi {x 2 )5y 2 {x2 1 xix 2 ). (3.6) 

Here <r(x 1 ,y 1 ;x 2 ,y 2 ) = (x 2 ,y 2 ;x 1 ,y 1 ). 

The Haar functional h : D(G) -» C and h* : D(G)* -» C are given by: 

h(f) = J dxf(x,e), h*(f) = J dyf(e,y). (3.7) 

Note that the Haar functional h* on D(G)* coincides with the co-unit e on D(G), 
and similarly h coincides with e*. Using these functionals, a positive-definite inner 
product on D(G) can be defined as: 



(/i,/2>= I dxdyh(x > y)f 2 (x,y) = h(h*f 2 *) = h*(f 1 *%). (3.8) 

To acquire a better understanding of the Drinfeld double, let us consider some 
of its sub-algebras. 



3.1. Sub-algebras of the Drinfeld double. The Drinfeld double is a twisted 
product of the algebra of functions on the group and the group itself. The subal- 
gebra A* of functions on the group described in the previous section is represented 
by the elements: 

f(g) -> 5 e (x)f(y) G D(G). (3.9) 

Then it is easy to check that the *-product coincides with the usual point-wise 
multiplication of functions: 

{S a fi*5 B f 2 )(x,y) = 6 e (x)f 1 (y)f a (v). (3.10) 

It is therefore clear that the *-algebraic structure on the algebra A* of functions 
on the group that we have described in section |2 is exactly the one obtained from 
the Drinfeld double dual algebra structure l|3.3|l when specialized to functions of 
the form S e (x)f(y). Note that the algebra A* with its pointwise multiplication can 
also be embedded into the Drinfeld algebra itself with its »-product. Indeed, as is 
easy to check, the "-product on functions of the form f(x)5 e (y) is just the usual 
pointwise multiplication. 

The group is represented by elements of the form: 

g ^ 5 g (y) E D(G). (3.11) 

The •-multiplication reduces to the usual group multiplication law: 



(An • 5 92){x,y) = S gig2 (y). 



(3.12) 
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3.2. Irreducible representations. Let us denote by Co, 9 6 [0, 7r] the conjugacy 
classes in SU(2), and by ge a representative of Cg that is in the Cartan subgroup 
U(\). The irreducible representations are labeled by pairs (6, k) of conjugacy classes 
and representations of the centralizer U(l). The carrier space V^ 8 '^ is: 

V {e > s) := {<j> : G -> C|0(arft.) = e - is «0(:r), Vx £ G,h = diag(e l « , e" i€ )} (3.13) 

and the action of an element / <G -D(G) is: 

(7T^ s \f)4>)(x) = J dz!{xg e x-\z)^x). (3.14) 

An orthonormal basis in V^' 8 ^ is given by the matrix elements: y / dim J T^ s (x), j > 
s, —j < m < j. We shall also use the bra-ket notation and denote the basis 
vectors by \jm). The carrier space V^ ' s ' can be decomposed into finite dimensional 
subspaces: 

y(M = @ j > s V^' s \ Vj e ' s) = Span{|jm)}. (3.15) 

3.3. Matrix elements and characters. Let us consider the matrix elements: 

(^• s) (f)YL' = (jm\^ e ' s \f)\j'm') = (3-16) 
ydlm~dim~7 J dxdz T? ns {y)f(yggy^ 1 1 z)T 3 m , s {z' 1 y). 

By definition, the matrix elements are in D(G)* . Using the pairing (-,-) we can 
identify them with functions on G x G. To this end, let us transform the above 
formula. Let us introduce, for each element x £ Cg of the conjugacy class Cg, an 
element B x £ G such that: 

B x g e B- x = x. (3.17) 
We will also need a 5-function 5g{x) picked on a conjugacy class 9, which is defined 



as: 



/ 



dgS e (g)f(g)= [ fixggx^dx. (3.18) 
Jgih 



IG/H 

Using these objects the matrix element can be written as: 



(^' S) (/))£>< - v/dhn.dmv jdxdyf(x,y)Sg(x)TUB,)Ti s (y- 1 B x ). (3.19) 



Therefore we define: 



(K (9 ' s) XL^y) ■■= V^J^fM^TUB^TCsiy-'Bz). (3.20) 

This is the main formula for matrix elements as functions on G x G. 
The character is obtained in the usual way as the trace: 

X {e > s) (x,y) : >>"'^; ,;„/.<•.//: -Wi^dhn/ei//,^ V B X ). (3.21) 

jm j 

To simplify this further let us use the fact that: 

J2^jTL(g) = e- is ^5(cos6). (3.22) 
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We have used the Euler parameterization here: g = g((j),9,tfj). Applied to (|3.21|1 
this formula implies that y = B x g^B~ x , or, in other words, that 13.21(1 contains the 
5-function S xy (yx). Thus, we have: 

X {6 ^\x,y) = 6e(x)S xy (yx) X s(B- 1 y- 1 B x ), (3.23) 

where Xs( n ) = s ls ^, n = diag(e^, e~ 1 ^). This is the character formula given in 16 . 

3.4. Decomposition of the identity. Now that we have obtained the characters, 
we can use them to decompose the identity element in the Drinfeld double. To this 
end, we will need the Weyl integration formula: 

f dgf(g)= [ -A(9) 2 f dx /(a^ 1 ). (3.24) 
Jg Jh/w n Jg/h 

Together with H3.18f) this implies: 

/ -A{9)H e (g) = l. (3-25) 
Jh/w n 

Let us now apply the summation over s to (|3.21|) . We get: 

f ^.A(9) 2 J2x ie ' s) (x,y) = Se(y). (3.26) 
Jh/w n ^-f 

4. Group field theory for the Drinfeld double 

From the algebraic formulation of Boulatov's theory in section [21 it was clear 
that the main object that is used in the construction of the model is the projector 
9 G A® 3 . Unlike the case of Boulatov theory, where there is essentially a unique such 
projector constructed from the identity element of A, in the Drinfeld double case 
there are several interesting "identity elements" that can be used in the construction 
of 9. We shall analyze several different possibilities. 

An explanation of why different possibilities can arise is as follows. The 9- 
projector that defines GFT should be constructed from a projector on D(SU(2)). 
Characters of irreducible representations are projectors. To construct a more gen- 
eral projector one can take various linear combinations of characters. Natural ex- 
amples are given by: (i) the sum of characters of all the representations of D(SU(2)), 
which gives the identity operator; (ii) the sum of characters of all simple represen- 
tations; (iii) the character of the trivial representation. One can more generally 
consider not characters but individual matrix elements, which are also projectors. 
In this paper we will consider the simplest case of the projector constructed from the 
identity operator, as well as another in a certain sense dual projector constructed 
from the element d e (x). 

4.1. Projector 9 constructed from the identity S e (y). Let us first consider a 
direct analog of Boulatov model. Thus, we construct a projector 9 from the identity 
element S e (y) on T>. To construct the projector we follow the same procedure that 
was used in section [2] Namely, let us first construct the propagator (or the kernel 
of the identity operator S e (y)): 



V* <g> V* 3 G = (K®id)A*l, 
G(x 1 ,y 1 ;x 2 ,y2) = S e (x 1 x 2 )6 e (y 1 y2 1 )- 



(4.1) 
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Let us define an operator h 3 : V* <g> D* <g> V* — > C via: 

M/1,/2,/3) :=h*(fi*h*f 3 )= (4.2) 
dz x dz 2 dz 3 dy fi{zi,y)f2{z2,y)h{z 3 ,y)5 e {z 1 Z2Z 3 ). 

The ^-projector is obtained by applying /13 to G® 3 : 

:= (id® 3 ® ft 3 ) G ® G <g> G. (4.3) 
Let us find an explicit expression for this projector. We have: 

9(xi,yi;x2,y 2 ;X3,y 3 ) = J dzxdz 2 dz 3 dy 5 e (x 1 z 1 )5 e (x 2 z 2 )6 e (x 3 z 3 ) x (4.4) 

8 e {z 1 z 2 z 3 )8 e {y^ 1 y)5 e {y^ 1 y)8 e {y^ 1 y) = J e (if x acJ 1 a;J 1 )^(yi, y 2 , 2/3), 

where 9(yi,y 2 ,y 3 ) is given by (|2.3U|) . An explicit verification shows that 9 is real 
9* = 9. We also have to check the projector property of 9. However, unlike the 
case of Boulatov model, the quantity (|4.4j) is not a projector. The "-product of two 
9 gives a divergent factor of rf 2 (|1.10(l that we have already encountered before 
in the discussion of the Ponzano-Regge model. We nevertheless proceed formally, 
and introduce a multiple of 9 as a projector. The following lemma is a statement 
to this effect. 

Lemma 4.1. The quantity i] 2 9, where ry~ 2 = S e (e) is a projector: 

( V 2 9) • { V 2 9) = V 2 9. (4.5) 

Proof. A proof is by verification. We have: 

(9»9)(x h y i ) = 6 e (xi 1 X2 1 x 3 - 1 ) (4.6) 
J d9 Y\. J dz i$g( z i) J dgS e (z^ 1 x^ 1 z 1 z 2 ' 1 x 2 1 z 2 z^ 1 x^ 1 z 1 )Y\_Sg(z^ 1 y i ) = 

i i 

(5 e (xi 1 x 2 1 x 3 ~ 1 )) 2 / dgdg'[[Sg(g~ 1 y l ) ^r/- 2 9(x t ,y t ). 

□ 

The "projector" so defined is the weakest point of our construction. A sceptic 
may argue that this quantity does not exist or is zero. We note, however, that 
both such a projector and the 77 factor would be perfectly well-defined had we been 
working with the Drinfeld double of the quantum SU g (2) at root of unity instead. 
Passing to the quantum SU g (2) is known to have the interpretation of adding a 
(positive) cosmological constant. Thus, as we see, unlike the Boulatov model itself, 
the model for the Drinfeld double of SU(2) is defined only formally. Later we shall 
consider another model which is well-defined. 

The quantity 9 can be expressed in terms of the i?-matrix. To state a result to 
this effect, we define the following graphical notation: 

R{4>® ip)(x 1 ,y 1 ; x 2 ,y 2 ) = (j)(xi,yi)il)(x^ 1 x 2 x u x^ 1 y 2 ) = II (4.7) 
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The opposite braiding is given by: 



The following result is confirmed by an explicit computation: 
Lemma 4.2. 



(4.8) 



(£ 'i 1 *...i n 1 M n ...iii,3; 1 x 1 1 . . . x n 1 xy) x (4.9) 



XiX^X yi) ■ 



Let us now take <j>(x,y) — S e (y), and take the Haar functional /i in the first 
channel. We get the following result: 

Lemma 4.3. 

Vi V 2 ¥„ 



-^..o/*n*<rw'*>. 



(4.10) 



We can use this result to give another description of the quantity 9: 
Lemma 4.4. 

1 1 1 



1 = 



(4.11) 



Let us consider the object (|4.1UI) . We will need two properties of such objects, 
which are referred to as the handleslide and killing: 

Lemma 4.5. A composition of two quantities H4.10|) satisfies the following "han- 
dleslide" property: 



Vi ¥ 2 ¥ k ¥ k+ ¥ k+2 V n 



Vi ¥ 2 ¥ k V k+ ,V k+2 ¥„ 



(4.12) 
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Proof. A direct computation of the left hand side gives: 

k 



JeOi 1 ■■■x k 1 ) / dgY[i>i(g 1 x i g,g 1 y i ) x 
i=i 

/n 
dg Yi Mg^xi^a^yi)- 



(4.13) 



The right-hand side is given by: 



S^x^ ...x n 1 )S e (x 1 1 ...x^) 

/k n 
dgdgY[A{g~ 1 9~ 1 x i gg,g~ 1 g' 1 yt) V'iO^iflSfl^V)- 



(4.14) 



i=l 



i=k+l 



A straightforward change of integration variable in the first part of the product in 
(14.1411 makes it identical to (|4.13|) and proves the lemma. □ 

Lemma 4.6. The following killing property holds: 



¥i ¥ 2 ¥„_, 



(4.15) 



Proof. We take ip n (x n , y n ) = S e (y n ), compute l|4.10|l . and apply the Haar functional 
h in the last channel. We get: 

« « ri — 1 ri — 1 

/ daJnJe^r 1 ---^ 1 ) / d 9 Y[^i(9~ 1 Xi9,9~ 1 yi) S e(g~ 1 ) = Y]_ i^fa, y t ) . (4.16) 



This proves the lemma. 



□ 



For completeness, let us also give a decomposition of the quantity 9 into charac- 
ters: 



Lemma 4.7. The following decomposition of the quantity 9 holds: 

9(xi,yi;x 2 ,y2;x 3 ,y 3 ) = / TT — - A(9i) 2 x 
JH/wfJ{ i" 

dim^ dim j2 dim j3 6>^ 2 2 ;/ 3 3 (xi , yi ; x 2 , y 2 ; x 3 , 2/3 ) , 



(4.17) 



where 



31,32 ,33 

^i(^i)^e 2 {x 2 )5e 3 {x 3 )5 e {x^ 1 x^ 1 x^ 1 )9 n ^ M {y 1 ,y 2 ,y 3 ), 



(4.18) 



where 9j i {yi) is given by <|2.31[1 . 

Proof. Let us use the character decomposition of the identity operator 5 e (y) : 

5 e (y)= [ -A^f^dim,^,*/), (4.19) 

JH/W n 
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where 

X^y) := £(^ s) )« m (x,y) = 6g(x) Xj (y). (4.20) 
We can now apply the formula l|4.10|l to a product of 3 characters, which results in 

Km . □ 

4.2. Projector constructed from 5 e (x). Another projector that we consider in 
this paper will be built from the operator l(x, y) = S e (x) £ D(SU(2)). This clement 
is not an identity: (1 • f)(x,y) — S e (x)h*(f). As we see, it plays the role of the 
element dual to the Haar measure on D(SU(2))*. We have put a hat above the 
symbol denoting this operator in order not to confuse it with the identity 5 e (y) in 
the Drinfeld double. It is clear that the ^-projector constructed from i according to 
the rules as described above is equal to ?7~ 2 1® 3 . Indeed, we take 3 operators 1 and 
enclose them with a strand with 1 inserted. The result (|4.9|) proves the assertion. 
It is obvious that a multiple of 6 constructed this way is (formally) a projector. 
Another important property of the operator 1 is that it satisfies the handleslide 
and killing properties. Let us state two lemmas to this effect. 

Lemma 4.8. When the projector P = 1 is inserted into the meridian link, the 
handleslide property holds. 

Proof. The left hand side is given by: 

k n 
X k+1 ■ ' ■ X n 

i ) n ( 4 - 21 ) 

i=l i=k+l 

The right-hand side is given by: 

k n 

^(zr 1 ...x^ l x^l l ...x~ 1 )'[[il) i (xi,yi)S e (x^l l ...x~ 1 ) Y[ AixiiVi)- (4-22) 

i=l i=k+l 

The handleslide property is obvious. □ 

Lemma 4.9. The quantity 1 satisfies the killing property provided the operator 
inserted into the longitude is a function of y only. 

Proof. A proof is by direct verification. □ 

Having discussed several different projectors that can be used in the construction 
we are ready to define the model. 

4.3. The model. There are now several possible models that can be formulated, 
depending on which ^-projector one uses. We shall proceed at the formal level, 
introducing factors of rj when necessary. Later the models we obtain will be placed 
on a more solid footing using the chain mail techniques. 

Let us consider a projector operator n 2 9. The basic field of the model is <f> £ 2?® 3 
that is real <fi* = 4>i an d one builds a projected field <j) via: 

4>:=n 2 6.0 = n 2 [7\f}~ (4.23) 
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In this graphical notation what is inserted on the "meridian" link and other 3 
strands depends on the model. In all cases rf'O is a projector. Using the graphical 
notation introduced, the action for the model is defined as follows: 




(4.24) 



To describe Feynman amplitudes generated by these models we shall introduce 
a notion of the chain mail. 



4.4. Roberts' chain mail. Using an idea of chain mail Roberts presented [T7] 
a very convenient description of the Turaev-Viro invariant. Here we remind the 
reader this notion, and show that the Feynman amplitude generated by our model 
is just a chain mail evaluation. 

Let us consider the Feynman diagram perturbative expansion of the model l|4.24f) . 
As the propagator of the model consists of 3 lines, each diagram is a collection of 
vertices (0-cells), edges (1-cells) and faces (2-cells) obtained by following each line 
till it closes. Each diagram D is an abstract one, that is not embedded in any space. 
Note that the data of a diagram D define a handlebody H(D), which is just the 
blow up of the graph of D. In other words, to obtain H(D) one takes solid balls 
(0-handlcs), one for every vertex of the diagram, and attaches them to each other 
by solid cylinders (1-handles), one cylinder for every edge. The handlebody H(D) 
is not embedded in any space. 

Now draw the curves defining the 2-faces on the boundary dH(D), and push 
them slightly inside of H{D). Let us then add the meridian curves Si for all the 
1-cells. The obtained collection of curves is called a chain mail link C(D) C H{D) 
of D. We now attach an appropriate (model-dependent) operator from T> to every 
component of the link C(D), and evaluate the link to obtain a value flC(D) G C. 
This evaluation procedure is as follows. One views the chain mail as a rule for taking 
a product of ^-projectors, one 9 for every 1-cell. The projectors are •-multiplied 
as specified by the 0-cells. Every closed loop corresponds to applying the Haar 
functional h. This evaluation rule will become more clear when we consider concrete 
examples below. We note that this evaluation is not an evaluation of the chain mail 
as embedded in some 3-manifold. It can be related to the more standard evaluation 
of the chain mail as embedded (as appears in Roberts' work |17|). see more remarks 
on this below. At the moment, however, there is no 3-manifold associated to any 
Feynman diagram. 

Proposition 4.10. For the model (|4.24|l the amplitude Z of a Feynman diagram 
D is equal to g v r] 2E QC(D), where V is the number of vertices and E is the number 
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of edges of D. Feynman amplitude of the model constructed using 1 is equal to 
g v nC(D). 

Proof. In the space of projected fields (j> = 9 • <f> the kinetic operator rfO is the 
identity operator. Its inverse is itself. Thus, the propagator of the model (|4.24(1 
is again r] 2 9, and the vertex is g times the operator depicted graphically in H4.24fl . 
This set of Feynman rules makes the first statement evident. The second case is 
similar, except that no factors of rj appear. □ 

4.5. Feynman diagrams and 3-manifolds. Each Feynman diagram D defines 
a 3D pseudo- manifold PM{D). It is easiest to obtain PM{D) by gluing together 
truncated tetrahedra an example of which is shown in Fig. 14.51 Thus, each vertex 
of D naturally corresponds to a truncated tetrahedron. Each edge of D defines a 
gluing of two truncated tetrahedra T,T': a large face of T is glued to a large face 
of T' . This way each (closed) Feynman diagram D defines a 3D pseudo- manifold 
PM(D) with boundary dPM (D). The pseudo- manifold PM(D) can be completed 
to a manifold M(D) if the boundary dPM(D) is a disjoint union of a number of 
spheres S 2 . In this case M(D) is formed by gluing to PM(D) a number of 3-balls. 
See ^H] for more details on this construction. 

Note, however that it may well happen for some of the Feynman diagrams that 
dPM(D) contains boundary components other than spheres. In this case one 
needs extra, structure to be added to the Feynman diagram to "convert" it into 
a 3-manifold. This extra structure is easiest to understand on the example of a 
toroidal boundary component in dPM(D). In this case, to convert the 3-manifold 
PM(D) with boundary into a closed 3-manifold one has to glue in a solid torus. 
However, there is no unique solid torus. A solid torus is specified by a closed curve 
on its boundary (torus) which is contractible inside the solid torus. Such a curve 
is in turn specified by two relatively prime numbers. This is exactly the extra data 
(for each toroidal boundary component) that are necessary to convert a Feynman 
diagram into a closed 3-manifold in cases when the manifold obtained by gluing 
in the two handles contains toroidal boundary components. When the boundary 
components are of higher genus, one needs even more data, which are a set of 3g — 3 
curves on the boundary {g is the genus of the boundary component in question) 
that are contractible inside the handlebody one glues in. These issues are delicate 
ones, and we would not like to divert the reader from the main theme by pursuing 
them any further. Let us just say that one either chooses (if possible) the symmetry 
properties of the model so that the boundary of PM(D) is always a collection of 
spheres, or, if the former is not possible, adds some extra structure to the model 
that tells one how to complete the Feynman diagrams into 3-manifolds. We will 
assume that one of the two is done, and proceed without worrying about these 
issues any more. 

Let us now explain a relation to the chain mail construction of the previous 
subsection is as follows. One obtains PM{D) by taking the handlebody H{D) and 
gluing to it a number of 2-handles (solid cylinders). The gluing is done as specified 
by the curves. Once PM(D) is obtained, one glues a number of 3-handles to ob- 
tain M(D). After this, using the "extra rules" discussed in the previous paragraph, 
one glues in a collection of 2-spheres or higher genus handlebodies, and completes 
PM{D) to a 3-manifold M{D) without boundary. The original Feynman diagram 
D is now embedded into M(D). The evaluation that we described above did not 
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need a representation of D as embedded into some 3-manifold, and in that sense 
was independent of any embedding. However, once the diagram is embedded into 
M(D), the evaluation that determines the Feynman amplitude can also be thought 
as the usual Reshetikhin-Turaev-Witten evaluation of D in M(D). The result of 
this evaluation does depend on the embedding of D in M(D). The amplitude one 
gets is thus that for a process described by the diagram D, happening in a particu- 
lar 3-manifold M(D). This should be contrasted with the usual situation in QFT, 
where the Feynman diagrams are not embedded in any space. The extra structure 
of faces and ordering of the faces that naturally follows from the group field the- 
ory expansion is exactly the structure that allows to think of D as embedded in a 
3-manifold. The 3-manifold that gets constructed from D with its extra structure 
is what describes the gravitational part of the degrees of freedom of the model. 

We now note that different Feynman graphs can lead to one and the same topo- 
logical manifold M(D). There are certain moves that relate different D that cor- 
respond to the same manifold M. As stated by the following theorem, when the 
operator that is inserted on the strands satisfies the handleslide and killing prop- 
erties, the chain mail evaluation is invariant under the moves and thus gives an 
invariant of M. 

Theorem 4.11. (Roberts) Assume that the operator O inserted on links forming 
C(D) satisfies the handleslide and killing properties. Assume that the trace of this 
operator is given by n~ 2 . The following multiple of the chain mail evaluation: 

CH(D) :=n 2no+2n3 nC(D), (4.25) 

where n and are the numbers of and 3-handles correspondingly, is an invariant 
of 3-manifolds m that CH(D) = CH(D') when M(D) ~ M(D') are 3-manifolds 
of the same topology. 

Proof. We shall only present a sketch, as a detailed proof is given in . Different 
Feynman diagrams that correspond to the same M are different handle decompo- 
sitions of M. A known theorem of topology states that different handle decom- 
positions can be related by a sequence of births of k, k — 1-handle pairs, and by 
handleslides. Invariance under handleslides is guaranteed by the the handleslide 
property of the operator inserted along the strands. To analyze the birth of 1-0 
and 3-2 handle pairs one uses the handleslide property, and is left with a single 
unknot with O inserted and not linked with the rest of the chain mail. This gives a 
factor of r}~ 2 absorbed by the prefactor of CH(D). A birth of the 2-1 handle pair 
is handled using the killing property. □ 
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Remark 4.12. Note that if the operator inserted along the meridians and longitudes 
satisfies the handleslidc and killing properties, and its trace is equal to unity, then 
a similar theorem holds, except that it is now the chain mail evaluation itself that 
is invariant. 

Corollary 4.13. The model constructed using the identity operator 5 e (y) gives 
invariants of 3-manifolds. 

Remark 4.14. There is another model that gives manifold invariants, namely the 
model constructed using 1. This model, however, requires a somewhat different set 
of formal prefactors. We shall describe it in details below. 

5. More general models 

In the previous section we have seen how group field theory for the Drinfeld 
double leads to the notion of chain mail in that Feynman diagrams of GFT are 
computed as the chain mail evaluation. Chain mail arises because operators that 
are used in the construction of the GFT model are projectors. One can therefore 
multiply the ^-projectors appropriately, and be left with a chain mail where there 
is just one meridian link per edge of a Feynman diagram D. Also importantly, 
the operator that is inserted in the strands is a projector, and after strands close 
combines to give a single operator in the longitude. 

One can consider a more general class of models not related to any GFT, but 
formulated directly in terms of the chain mail. Thus, one takes a chain mail that 
corresponds to a graph D, inserts one species of operators into the meridians, some 
other species into the longitudes, and evaluates the resulting link. The operators 
inserted do not have to be projectors anymore, and this is what makes such models 
different from the GFT ones. However, as we shall see, these more general models 
admit a physical interpretation and are of interest. We could have directly started 
from the notion of chain mail and formulated all models correspondingly. However, 
we believe that the GFT description we have given is essential in that it clear 
shows what is and what is not possible in the GFT framework. The GFT models 
described are also of importance in view of possible generalization to algebras other 
than D(SU(2)). Thus, for instance, the analog of Boulatov model for the quantum 
group SU 9 (2) has not been formulated. It is straightforward to do so using our 
algebraic formulation described above. 

5.1. Formulation of the models. Let us define a set of models as follows. Con- 
sider a chain mail evaluation in which an operator P = 5 e (x)P(y) is inserted in 
all the longitudes. The function P(y) that appears as part of this operator is not 
required to be a projector. Instead, we will just require P(y) to be normalized so 
that J dyP(y) = 1. Let us consider the following model. 

Definition 5.1. The amplitude of the model is obtained by inserting the identity 
S e (y) into the meridian links, and 5 e (x)P(y) into the longitudes. The amplitude is 
a (formal) multiple of the chain mail evaluation: 

CH'(D) = Tf 2ni QC(D). (5.1) 

Here denotes the evaluation, and the factor of i] 2 to the power of the number of 
1-handlcs is introduced for reasons that will become clear below. 
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Let us prove some properties of the objects that the both models are constructed 
from. First let us consider a number of strands enclosed by the P operator. Such 
an object is given by: 

P(x^ 1 ...x- 1 )]lrl> i (xi,y i ). (5.2) 

i 

Importantly, there is no longer a handleslide property 14.121 for objects (|5.2|l . as is 
easy to verify. However, the killing property still holds. Indeed, we take tp„ = 6 e (y), 
and apply the Haar functional in the last channel. We get: 

/n— 1 
dx n P(x^ 1 . . . x^ 1 ) Y[ i>i{xi,yi)5 e (e) = (5.3) 
»=i 

/n— 1 n— 1 

dgP(g) Y[ $i{ x iiVi) = II fafa'Vi)' 

i=l i=l 

where we made a change of integration variable to arrive to the second expression. 
Using these facts, it is easy to prove the following assertion. 

Proposition 5.2. Amplitudes given by (15.1(1 are invariant under 0-1 and 2-1 han- 
dle pair births /deaths, as well as under 1-handle slides. 

Proof. Proof is same as that of Roberts' theorem of the previous section. The 
factor of r] 2ni in (|5.1|l is necessary to guarantee the 0-1 handle pair births/deaths 
invariance as well as the invariance under the 2-1 handle pair births/deaths. Note 
that in the case considered by Roberts, see l|4.25() above, one needs two prefactors 
to guarantee the 0-1 and 2-3 moves. In our case we do not require the 2-3 moves, 
but have an additional factor of rf 2 appearing in the 1-2 moves. This makes the 
correct prefactor to be rj 2 to the power of ri\ not no as in case of l(4.25|l . □ 

5.2. A model giving 3-manifold invariants. Let us note that the case P(y) = 1 
is special, because in this case there is the 2-handleslide property. In this case the 
chain mail evaluation gives 3-manifold invariants. We have the following (almost 
trivial) corollary. 

Corollary 5.3. Evaluate the chain mail by inserting the identity operator S e (y) in 
all meridian links, and the operator 5 e {x) into all longitudes. The quantity (|5.1H is 
a topological invariant in that it only depends on the topology of M(D). Moreover, 
its value is one for any manifold M(D). 

5.3. Interpretation of 2-handleslide invariance absence. The example of the 
previous subsection shows that there are models that are interesting but not in 
the class of GFT models considered in the previous section. It is not hard to give 
a more flexible definition of GFT that would cover more examples. We shall not 
attempt this however, concentrating instead on the chain mail definition from now 
on. 

General models defined via chain mail do not produce 3-manifold invariants 
because there is no handleslide property for the operator inserted in the longitudes. 
Still, as we shall presently see, from the point of view of point particle theory these 
models are quite natural. To see this, let us analyze what happens if there is no 
handleslide property. Meridian links correspond to 1-handles, and the longitudes 
correspond to 2-handles. There is the handleslide property for the identity operator, 
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see lemma l4*.12l Thus, the chain mail evaluation Q,CH(D) is invariant under the 1-0 
handle pair births and deaths. The killing property guarantees an invariance under 
the 2-1 pair births and deaths. There is also invariance under the handleslides along 
1-handles. Thus, non-invariance comes only from two sources: there is no invariance 
under 2-handle slides, and there is no 3-2 handle pair birth-death invariance. Now 
recall that each D is the dual skeleton of some triangulated 3-manifold M(D). As 
we have discussed in the introduction, we would like to interpret the triangulation 
as a Feynman diagram as well. Edges of this Feynman diagram are in one-to-one 
correspondence with 2-handles, and vertices with 3-handles of M(D). Absence of 
invariance under moves involving 2-handles just means that the amplitude we get 
is not invariant under changes in the Feynman diagram. The absence of invariance 
under 3-2 pair birth-death means that the Feynman amplitude is not invariant under 
a vertex being replaced by two vertices connected by a new edge. Similarly, the 
absence of the 2-handleslide property means that one cannot move edges through 
vertices along other edges. All these moves change the Feynman diagram, and the 
fact that there is no invariance is very natural from the point of view of point particle 
field theory. Thus, in this case the model leads not to 3-manifold invariants, but 
to Feynman amplitudes of the dual point particle field theory. To get most general 
point particle models we will have to work with theories of this class. 

6. Evaluation and interpretation of the models 

We have considered a set of models. Some of them were shown to give invariants 
of 3-manifolds, some other only give amplitudes for dual Feynman graphs. To 
analyze the structure of the amplitudes arising, and to give it a geometrical and 
point particle interpretation, let us use the kernel representation. 

6.1. Kernel representation. Recall that in the kernel representation one asso- 
ciates an operator f(x,y) € T> its kernel defined as (k ® id)A*/. A simple compu- 
tation gives: 

f(x, y\x, y) = f(y~ 1 x~ 1 y, y~ 1 y)6 e (xx). (6.1) 

The kernels should be multiplied using the ^-product and the h* Haar functional 
should be taken. In the case of algebra of functions on the group this operation is 
just a point-wise multiplication of the propagators with an integral taken. In the 
case of the Drinfeld double the structure is more involved. To display it, we note 
that: 

h*(fi*f 2 ) = J dxdy f 1 (x,y)f 2 (x-\y). (6.2) 

Thus, in order to compose kernels of two operators, one should take one kernel and 
multiply it by the other kernel with x replaced by x~ 1 , and then integrate over x, y. 
Equivalently, one can replace x by x" 1 in (|6.1I) . It is this quantity that we shall 
refer to as the operator kernel: 

Kf(x, y\x, y) := f(y~ 1 xy, y~ 1 y)6 x (x), (6.3) 

One can now multiply the kernels Kf in the usual way, and integrate over the gluing 
variables. Thus, it is easy to check that: 

J dpdq K fl (x,y\p,q)Kf 3 (p,q\x,y) = K^.f 2 {x,y\x,y). (6.4) 
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However, the kernel formalism leads to a singularity for the Haar functional h. 
Indeed, we have: 

J dxdy K f (x, y\x, y) = S e {e)h(f) = v^Hf). (6.5) 

Thus, in order to reproduce the correct Haar functional the trace of the kernel 
should be renormalizcd by the already familiar quantity rj 2 (|1.10|l . 

6.2. Group field theory Feynman diagram analysis. Let us first apply the 
kernel technique to the model that was constructed from the identity operator 
S e (y) and that gives 3-manifold invariants. In this subsection we group quantities 
according to GFT Feynman graphs D. In the next subsection we perform a "dual" 
analysis, in which a point particle interpretation is more clear. 

It is clear that the amplitudes can be computed as follows: one should take 
the kernel for each vertex of D, and compose them together as specified by the 
diagram. This is due to the fact that the propagator of the model is proportional 
to 9, and the vertex in (|4.24(l is invariant under multiplication by 9 from any of 
the 4-possible directions, see Q4.24[l . Let us analyze this vertex kernel and give its 
geometrical interpretation. 

First we need an expression for the ^-projector kernel. Applying l|6.3|l to (|4.4() 
we get: 

Kg(xi, J/i ; X 2) 2/2, ^3, V3\Xl, yV,X2, 2/25 X 3 , 2/3) = $xi (xi)S X2 (x 2 )S X3 (x 3 ) (6.6) 

Se(yT 1 xmy2 1 x2y2y3 1 x 3 y 3 )9{yi 1 y 1 , y^m-, y^m)- 

Let us now find the vertex kernel. As in section [2 we shall enumerate 4 channels 
of the vertex by indices i,j = 1, ... ,4. The vertex kernel is then a function of 
variables Xij,yij, where Xij ^ %ji,yij 7^ 2/j'i- The kernel is obtained by taking a 
composition of 4 kernels Hfi.6|l . We have: 




where x%j = Xji, jjij — 2/ji- It is straightforward to do the x integrations. Each kernel 
H6.6|) contains a 8 x {x) factor, which makes the result proportional to S Xij (xji). Thus, 
the vertex kernel is actually a function of 6 variables x\j = Xj^. 

V(x ij ,y ij )= ! ! 5 e (yr^Xi j y ij yT^x ik y ih yr l 1 x a y a ) (6.8) 

Here jjij = jjji, so the integration is taken over 6 variables. We have also introduced 
a notation: k, I) for a quadruple of integers: i j ^ k ^= I. The last y integral 
in (|6.8|l can be taken at the expense of introducing another 4 variables <?j. Indeed, 
let us use the formula H2.30fl for each of the 4 #-functionals in (|6.8|) . We get: 

V(xij,yij) = Y[ 8 e {y~ j l x lJ y lJ yl k 1 x lk y lk yl l 1 x a y l i) J JJdffi 5 e {y lJ g l gJ 1 yj i 1 ). (6.9) 

(i,j,k,l) J i<j 

We note that the last term in this expression is just the vertex kernel of the Boulatov 
model. Thus, the modification of the vertex as compared to the Boulatov model 
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case is in an appearance of the 4 additional <5-function terms, and in the fact that 
the kernel depends on 6 variables Xij . 

One can re-interpret this structure by saying that the only modification as com- 
pared to the Boulatov model case is that the edge contribution changed: there is 
now an extra 5-function for each edge of D. Let us change the notation to display 
the structure arising more clearly. Let us denote the vertices, edges and faces of D 
by tilded letters: v, e, / correspondingly. We leave the untilded letters for elements 
of the triangulation T to which D is dual. Let us now denote the y-variables as 
y~ f, where e is an edge, and / is a face of D (each edge e is shared by exactly 3 
faces). Let us denote the Xij variables by xi. Indeed, because there is a (5-function 
imposing x t j — Xji, all x- variables belonging to edges around a face are equal, so 
there is just one x-variable per face. Thus, the edge factor becomes: 

'■^l !>..)■>■ j' !>..<> (6-10) 

/ 

for each edge e of D. The vertex, on the other hand, is given by the same expression 
as in the case of Boulatov model: 

J II dff5 < e ~ S e(ye,f9v,e9^ly~j), (6.11) 

where for each pair e, e! of edges that emanate from vertex v f is the face that shares 
both e, e'. Note that the (/-variables are different at different vertices g$,e ^ 55', e- 
To obtain the amplitude one multiplies the edge and vertex kernels, and integrates 
over the y~ ?, Xj variables. 

Already at this stage an interesting geometric interpretation of all the variables 
is possible. Indeed, we have already encountered the g% t g variables before. Recall 
that each Feynman diagram D is dual to a triangulated 3-manifold, and therefore 
each vertex v is dual to a tetrahedron t. Thus, all the variables in l|6.11f> are those 
of a single tetrahedron. Let us denote the triangulation to which D is the dual 
graph by T. We denote the elements of T by untilded letters. Thus, v, e, /, t are 
vertices, edges, faces and tetrahedra of T correspondingly. We see that g^ : e are 
variables of a tetrahedron t dual to v and correspond to faces / of T . Each g^g has 
the interpretation of the holonomy of the connection across the face /. There is 
a similar variable in the neighboring tetrahedron t' and the total holonomy across 
the face / is given by the product: g/ = (g t j)~ l g t f . As we have already discussed 
in the introduction, the product (g t j : )~ 1 g t j:, describes a rotation of one face / into 
the other /' and thus carries information about the dihedral angle between these 
faces. 

The variables y~ j are more interesting and this is the first time that we have 
encountered them. Using the duality, we can also write them as y e j- There is 
12 of them for each tetrahedron. It is natural to interpret them in terms of a 
certain other truncated tetrahedron. Indeed, let us introduce the dual tetrahedron, 
and cut off its vertices by surfaces parallel to the faces of t to obtain a truncated 
tetrahedron t. A compound object made of t and i is shown in Fig. 16.21 The 
truncated tetrahedron t has 12 vertices, and we shall interpret 12 variables y e .f as 
describing "positions" of these vertices. The meaning of each (5-function in (|6.11(l 
is then as follows: it imposes the constraint that the rotation of a vertex y e j into 
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Figure 2. A compound of a tetrahedron t with a truncated tetra- 
hedron i. 

y e ,f around the edge e is the same as rotation of the face / into /' as described by 
g- variables. 

Thus, we only have the variables Xj uninterpreted. By duality, they are variables 
x e that have to do with the edges of t. Tetrahedron t is truncated, and we interpret 
x e as describing rotations of the small faces of the truncated tetrahedron t into 
one another. Thus, interpretation of all of the variables in (|6.9|l is in terms of 
geometrical quantities associated with two truncated tetrahedra t, t inscribed into 
one another. 

For the convenience of the reader, let us re-write the above edge and vertex 
factors using the triangulation T notations. We have: 

MlI^KV/) (6-12) 

e 

for each face / of T. Here the product in the argument of the (5-function is taken 
over 3 edges that form the boundary of face /. The tetrahedron factor is given by: 

/rWII S e {y e jg}{g),)- l y-},). (6.13) 
/ fj'et 

The interpretation that we would like to propose for the total tetrahedron am- 
plitude (|6.9|l is that of a wave function of a pair of dual to each other truncated 
tetrahedra. To obtain an amplitude for a manifold one multiplies these wave func- 
tions and integrates over the gluing variables x, y. We would like to note that 
there is a similarity in the structure of the vertex found and the expression for the 
6j-symbol obtained in |20| . In both cases the quantity in question has to do with 
a truncated tetrahedron, and is constructed using propagators for long edges (the 
second set of 5-functions in (|6.9|) ). and a set of factors for the small faces (the first 
set of 5-functions in Hfi.Hjl ). Thus, the structure we have found is probably pertinent 
to a very general class of models. 

If not for the edge factors (|6.12l) , one could easily integrate over the y variables 
with the result being a product of ^-functions, one for every edge e of T, and re- 
quiring that the holonomy around e is trivial. With factors (|6.12() the integration is 
not straightforward. The total amplitude (|6.9|1 is one where tetrahedron constraints 
discussed in the introduction are taken into account. 

6.3. Dual graph analysis. Above we have represented the chain mail evaluation 
ilCH(D) as a composition of vertex kernels (|6.9(l . or the edge H6.1()(l and vertex 
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l|6.11|l factors. Edges and vertices are those of D. We would now like to group the 
factors according not to D but to a triangulation T whose dual graph is D. 

Thus, Feynman diagrams D have an interpretation of the dual skeleton of some 
triangulated 3-manifold. Let us denote this triangulation by T. Let us consider 
vertices v of this triangulation. The edges e of T are formed by parts of the meridian 
links of CH(D). The number of strands forming each edge of T is not fixed, and 
is equal to the number of faces of T sharing this edge. Strands forming the edges 
of Feynman diagram D, and whose closure defines the faces of D, now enclose the 
strands forming the edges of T. Thus, we have a very similar structure to that 
considered above, except that it is a dual one, and that the number of strands 
enclosed by a link is now not fixed. To evaluate CH(D) we have to insert the 
identity operator in each strand. Thus, we insert in each edge e of T the following 
operator: 

S e (x^ ...x- 1 ) f dxUSeix^yi). (6.14) 

As before, it is convenient to introduce the kernel (|6.3() : 

^(yr^rV • ■ ■Vn lx n 1 yn) j dx '[[S x (y^ 1 y t )S Xt (x l ). (6.15) 

As in the previous subsection we can take a composition of such kernels around 
every vertex of T by integrating over the i, y variables. Taking into account the 
$xi{%i) functions, we see that there remains a single variable Xf for every face. 
The y- variables can be denoted by y e j- Indeed, there is one such variable for each 
strand forming each edge of T, or, equivalently, for each face / sharing edge e. This 
vertex kernel is found to be: 

II ^(Q VeJ^Vej) II 6 e(yejXv,eX~^ye,j)- (6.16) 
e / e,e' 

Here the second product is over pairs of edges that emanate from v, and / in 
the argument of the second set of 5- functions is the face shared by both e, e'. 
Equivalently, one can say that second first product is over the faces that touch 
v. To obtain the chain mail evaluation one should multiply the kernels (I6.16|) 
and integrate over all the variables x v ^ e ,Xf,y e j. Note that the vertex Ij6.16|) is 
essentially the same as ()6.9JI . except that the dual quantities are used. 

Let us discuss a geometrical interpretation of the quantities involved in Ij6.16|l . 
The quantity Xf is interpreted as the holonomy across the face /. The quantity 
y e j encodes the dihedral angles of a tetrahedron. Namely, 2/ e ,/yj}< describes the 
rotation of face / into /' around edge e. We see that arguments of the first set of 
(5-functions in (|6.16|) are just the product of holonomies across the faces that share 
e with the group elements that describe rotations of / to the next face /'. Thus, 
the argument of the <5-function for edge e is the total rotation that one gets by 
going around e. This means that there is no deficit angle allowed! Thus, in spite of 
the fact that the model was constructed using D(SU(2)), not just SU(2), the model 
does not seem to describe point particles. We shall return to this conundrum below. 

For now let us note that an interpretation of the second set of (5-functions in 
(I6.16J1 is also possible. Indeed, we see that the quantity x v ^x~\, describes the 
rotation of edge e into e', with the vertex v as the center of rotation. Each S- 
function guarantees that y~jy e 'j is the same as the rotation that sends e to e! . 
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6.4. General models and their point particle interpretation. As we have 
just seen, the topological GFT model constructed from the identity operator S e (y) 
on the Drinfeld double does not describe point particles. Earlier we have introduced 
a more general set of models in which a function 8 e {x)P(y) is inserted into the 
merideans. Now we would like to show that these models do describe point particles 
and P(y) is the particle's propagator in the momentum space. 

We have seen in the previous subsection, the point particle interpretation, if 
any, is most visible in a description in which one groups all the projectors around 
vertices of T. Let us give this description. As before, each edge of T is formed by 
a number of strands with the identity operator inserted. However, the operator on 
the link that encloses the edge is now P. Therefore, each edge of T corresponds to 
the following operator: 

p(*r 1 ...onw- (6 - 17) 

i 

These operators now have to be multiplied (using the ^-product) as specified by the 
fat structure, and the Haar functional has to be applied in each strand to get the 
amplitude. It is easiest to understand what the result is by noticing that each of 
these edge operators is actually an element of the ideal (A*)® n . Indeed, functions of 
the form f(x)5 e (y) form an ideal in D(G), isomorphic to the algebra A* with point- 
wise multiplication. Therefore, the result of the ^-product of all these operators 
is the point-wise product of the functions P{x^ . . -x' 1 ) times the product of the 
5- functions 5 e (y). It is now clear what the evaluation is. Let us first discuss what 
happens with the y-dependence. Applying the Haar functional h in each strand we 
get an infinity 8 e (e). These infinities are cancelled by the prefactor rf ni that we 
introduced in (|5.1|l . One could, actually, avoid this infinities altogether by a slight 
modification of the evaluation procedure. Thus, since the result in each strand is 
always in A* ideal of D(G), let us instead of applying the Haar functional of D(G) 
apply the one of A* . The result of this evaluation procedure is finite (for sufficiently 
nice propagators for which the group integrals would converge). 
The result of this evaluation procedure is given by 

j f 

Here the first product is over the faces of the diagram, and Xf is a group variable 
that gets associated to each face. The second product is over the edges of the 
diagram, and the product of a;- variables in the argument of the propagator includes 
these variables for the faces /i , . . . , /„ shared by that edge. It is clear now that the 
function P{y) plays the role of the particle propagator, for the above answer for 
the amplitude is as expected, see e.g. PE]. 

Thus, we have arrived at the expression we started from in the Introduction. 
Indeed, the amplitude is just the Feynman amplitude in the momentum represen- 
tation, with the propagator evaluated at the product of face variables around an 
edge. The momentum conservation constraints are solved by introducing the face 
variables. The tetrahedron constraints that we discussed in the Introduction can be 
omitted because the integrand is invariant under the action of the gauge transfor- 
mations at vertices of the dual triangulation, as can be easily checked. Moreover, as 
the proposition 15 . 21 states, the Feynman amplitude is to a large extent independent 
of the extra "fat" structure that is added to the diagram in order to define it. 
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There are two interesting "limiting cases" to consider. The case P(y) = S e (y) 
corresponds to having no particles at all. This gives the topological Ponzano-Regge 
model, if one introduces an additional prefactor rj 2n ' i that is required to guarantee 
the invariance under the 2-3 moves. The other case is P(y) = 1, which gives 
another (trivial) topological model, the one that gives one for every manifold as the 
amplitude. However, one can decompose the function P(y) = 1 into (5-functions 
on conjugacy classes as in l)3.25[l . This shows that even this trivial model can be 
thought of as containing particles of all possible masses in it, but becoming trivial 
after the integral over the particles' mass is taken. 

7. Discussion 

Let us recap the main points of this work. We have considered the group field 
theory for the Drinfeld double and constructed an exact analog of the Boulatov 
model. The Feynman expansion of the model leads to triangulated 3-manifolds 
weighted with a manifold invariant. The invariant is the "Ponzano-Regge" model 
for the Drinfeld double. By reducing this Feynman amplitude to the chain mail 
evaluation we proved that this Drinfeld double Ponzano-Regge amplitude is indeed 
an invariant of 3-manifolds (when dressed by an appropriate power of the formal 
quantity rj 2 ). We have also found that this model does not describe point particles. 
Indeed, one of the factors appearing in the amplitude for this model is a (5-function 
for the holonomy around each triangulation edge. This is a bit disappointing, be- 
cause our original motivation for considering the Drinfeld double GFT was precisely 
to obtain a model containing point particles. Thus, in a certain sense, our original 
motivation has failed. 

Nevertheless, motivated by the GFT considerations, and especially by the chain 
mail construction, we have introduced a more general set of models that are defined 
using not GFT but directly in terms of the chain mail. We have considered only 
one such model, but others are possible. The data that goes into this model is an 
assignment of a propagator (function on the SU(2) group manifold) to each edge of 
the Feynman diagram. This propagator can be chosen different for different edges, 
as long as it is appropriately normalized. The resulting Feynman amplitude is 
invariant under moves that do not change the Feynman graph and the 3-manifold 
topology. This model describes point particles. In fact, it is exactly the naive 
model envisaged already in the Introduction. Our results on its invariance under 
the topological moves show that this model is well-defined and to a large extent 
independent on the "fat" structure that was introduced to define it. The model 
models become topological if one choses P(y) = 1. In this case it seems trivial as 
the amplitude it gives to any manifold is one. At the same time, even in this case 
the model does describe point particles. One just has to interpret the operator 1 
as the integral over particle's mass, or using physics terminology, the integral over 
the internal momenta of particles. 

The model we have considered (the one constructed with a general propagator 
P inserted into the longitudes) has an interpretation as describing pure 3d gravity 
coupled to point particles. As we have discussed in the introduction, there is 
another possible interpretation. Namely, one can view this model as describing 
point particles only, but with the backreaction on the metric taken into account. We 
have not discussed what would be the corresponding theory in the position space. 
Let us only note that non-commutativity of momentum implies non-commutativity 
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in the position space as well. It would be of interest to formulate point particle 
theories described here in the position space, and understand what kind of action 
principle is responsible for them. 

Now that one has a model for 3d quantum gravity with point particles, various 
quantities of physical interest can be computed. However, certain additional de- 
velopments are necessary in order to achieve this. First it is important to be able 
to compute not only closed but also open Feynman graph amplitudes. This would 
give correlation functions. Second and more important, one should find an analog 
of LSZ reduction formula of field theory that expresses scattering amplitudes in 
terms of field correlation functions. We shall not attempt this in the present paper. 

A remark is in order about the topological invariant defined by our particle 
model. As we have already mentioned, this is Z — 1 for any manifold. However, 
the amplitude for each Feynman diagram (with a fixed set of conjugacy classes 9 e 
on the edges) is non trivial. Moreover, an amplitude for an open Feynman graph, 
which describes particle scattering and which we have not considered in this paper 
is also not trivial, even when one integrates over the internal momenta. This shows 
that, in spite of seeming triviality, this model is physically interesting. 

Let us now return to the topologically invariant model constructed from the 
identity operator S e (y). Several comments are in order. First, it would be inter- 
esting to understand how different this model is from that of Ponzano-Regge. The 
original Ponzano-Regge model imposes the constraints that the holonomy either 
around edges or dual edges is trivial. The model constructed in this paper im- 
poses both sets of constraints at the same time. For the case of the double of a 
quantum SU 9 (2), which makes the model well-defined (by rendering r/~ 2 finite), 
one expects the Drinfeld double "Ponzano-Regge" model to be equal to the usual 
Ponzano-Regge model amplitude squared. It would be interesting to understand 
the interpretation of the classical group Drinfeld double case better. 

Another remark is that, as is clear from our construction, the gauge symmetry 
8<f> = ip, • tp = of the group field theory is extremely important. From the chain 
mail description it is clear that topological invariance of the model is intimately 
related to the fact that ^-projector is used as the propagator of the model. On 
the other hand, we know that topological invariance of 3d gravity is due to large 
amount of gauge symmetry present in it. More precisely, diffeomorphisms are so 
strong in 3d, that they render the theory topological. All this strongly suggests 
that there is a link between the group field theory and gravity gauge symmetries. 
In the context of 2d gravity matrix models the relation between gravity and a 
matrix model is an instance of open-closed string duality, in which one and the 
same Feynman amplitude has two different interpretations. We expect that the 
relation between group field theory and 3d gravity explored here is of a similar 
nature, and that the group field theory gauge symmetry becomes diffeomorphism 
symmetry on the gravity side. It would be of considerable interest to explore all 
this in more detail. 

It is clear that the GFT construction that we have presented is very general, and 
can be applied to algebras other than the Drinfeld double. The Drinfeld double of 
SU(2) is described ^BJ as the twisted product of the algebra of functions on SU(2) 
with its dual, which is also algebra of functions on SU(2). The Drinfeld twisted 
product construction can be applied to groups other than SU(2). Another im- 
portant for physical applications example is given by the quantum Lorentz group, 
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which is just the twisted product of the quantum SU(2) with q real and its Pontrya- 
gin dual, see ^Hj- Our construction can be applied to this case. We expect that in 
this case there is not one, but several possible models that lead to 3-manifold invari- 
ants, depending on which "identity" operator is used. These models should have 
the physical interpretation as describing point particles coupled to quantum gravity 
with either positive cosmological constant and Lorentzian signature or with nega- 
tive cosmological constant and Euclidean signature. It would be of great interest 
to find all these models. 
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Appendix: Some useful formulas 

In the main text we often use the integral of two and three matrix elements. For 
two matrix elements we have: 

dg Tl n (g)T r C n ,(g) = Jj^wW (7.1) 

In graphical notation this takes the following form: 

d 9 fi x ; (I) = ' ( 7 - 2 ) 









r 


J dim 






and 

/ dg (g) (g) (g) = | | | (7.3) 
The intertwiner that is used in the last formula is chosen to be normalized so that: 



= 1 (7-4) 



whenever the 3 spins satisfy the triangular inequalities, and zero otherwise. 
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